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Algebraic Combinatorics of Magic Squares 

Abstract 

The problem of constructing magic squares is of classical interest and the first known magic square 
was constructed around 2700 B.C. in China. Enumerating magic squares is a relatively new problem. 
In 1906, Macmahon enumerated magic squares of order 3. 

In this thesis, we describe how to construct and enumerate magic squares as lattice points inside 
polyhedral cones using techniques from Algebraic Combinatorics. The main tools of our methods 
are the Hilbert Poincare series to enumerate lattice points and the Hilbert bases to generate lattice 
points. With these techniques, we derive formulas for the number of magic squares of order 4. We 
extend Halleck's work on Pandiagonal squares and Bona's work on magic cubes and provide formulas 
for 5 X 5 pandiagonal squares and 3x3x3 magic cubes and semi-magic cubes. 

Benjamin Franklin constructed three famous squares which have several interesting properties. 
Many people have tried to understand the method Franklin used to construct his squares (called 
Franklin squares) and many theories have been developed along these lines. Our method is a new 
method to construct not only the three famous squares but all other Franklin squares. We provide 
formulas for counting the number of Franklin squares and also describe several symmetries of Franklin 
squares. 

Magic labelings of graphs are studied in great detail by Stanley and Stewart. Our methods 
enable us to construct and enumerate magic labelings of graphs as well. We explore further the 
correspondence of magic labelings of graphs and symmetric magic squares. We define polytopes of 
magic labelings of graphs and digraphs, and give a description of the faces of the Birkhoff polytope 
as polytopes of magic labelings of digraphs. 
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Foreword 



This thesis is a story of fusion of classical mathematics and modern computational 
mathematics. We redefine the methods of solving classical problems like constructing 
a magic square and exhibit the power of computational Algebra. For example, only 
three Franklin squares were known since such squares were difficult to construct. 
Our results make it possible to construct any number of Franklin squares easily. A 
detailed description of our methods is given in Section 11.21 More significantly, our 
study leads to an elegant description of the faces of the Birkhoff polytope as polytopes 
of magic labelings of digraphs. The main theorems are introduced in Chapter ^ and 
the subsequent chapters explore the details of these results. 
This thesis is based on my three papers 

• Polyhedral cones of magic cubes and squares (joint work with J. De Loera and 
R. Hemmecke) 

• How many squares are there, Mr. Franklin?: Constructing and Enumerating 
Franklin Squares jS]; 

• Magic graphs and the faces of the Birkhoff polytope jl]. 
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Chapter 1 

Introduction 

1.1 Magic Squares. 

Miracles seem very wonderful to the people who witness them, and very simple to the people 
who perform them. That does not matter: if they confirm or create faith they are true 
miracles. 

- George Bernard Shaw. 

A magic square is a square matrix whose entries are nonnegative integers, such that 
the sum of the numbers in every row, in every column, and in each diagonal is the 
same number called the magic sum. The study of magic squares probably dates back 
to prehistoric times 0. The Loh-Shu magic square is the oldest known magic square 
and its invention is attributed to Fuh-Hi (2858-2738 B.C.), the mythical founder of 
the Chinese civilization [7j. The odd numbers are expressed by white dots, i.e., yang 
symbols, the emblem of heaven, while the even numbers are in black dots, i.e., yin 
symbols, the emblem of earth (see Figures fl . II and fT^ A) . 

Like chess and many of the problems founded on the figure of the chess-board, the 
problem of constructing a magic square probably traces its origin to India A 
4x4 magic square is found in a Jaina inscription of the twelfth or thirteenth century 
in the city of Khajuraho, India (see Figure E2lB). This square shows an advanced 
knowledge of magic squares because all the pandiagonals also add to the common 
magic sum (see Figure IT^ . 
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Figure 1.1: The Loh-Shu magic square [J]- 
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Figure 1.2: Loh-Shu (A), Jaina (B), and the Diirer (C) Magic squares. 



From India, the problem found its way among the Arabs, and by them it was 
brought to the Roman Orient. It is recorded that as early as the ninth century magic 
squares were used by Arabian astrologers in their calculations of horoscopes etc. 
which might be the reason such squares are called "magic" Their introduction 
into Europe appears to have been due to Moschopulus, who lived in Constantinople 
in the early part of the fifteenth century. The famous Cornelius Agrippa (1486-1535) 
constructed magic squares of the orders 3, 4, 5, 6, 7, 8, 9, which were associated 
by him with the seven astrological "planets"; namely, Saturn, Jupiter, Mars, the 
Sun, Venus, Mercury, and the Moon. A magic square engraved on a silver plate was 
sometimes prescribed as a charm against the plague and a magic square appears in 
a well known 1514 engraving by Albrecht Diirer entitled Melancholia (see Figures 
11.21 C and II. 3p . Magic squares, in general, were considered as mystical objects with 
the power to ward off evil and bring good fortune. The mathematical theory of the 
construction of these squares was taken up in France only in the seventeenth century, 
and since then it has remained a favorite topic of study throughout the mathematical 
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Figure 1.3: Magic square in a 1514 engraving by Albrecht Diirer entitled Melancholia 

world. See |Zj, 05 111], or to read more about the history of magic squares. 

Constructing and enumerating magic squares are the two fundamental problems 
in the topic of magic squares. Let M„(s) denote the number of n x n magic squares 
of magic sum s. In 1906, MacMahon enumerated magic squares of order 3: 

+ |s + 1 if 3 divides s, 

9 3 

otherwise. 

Later, MacMahon Anand et al. jB], and Stanley jB?!, considered the problem 
of enumerating semi-magic squares which are squares such that only the row and 
column sums add to a magic sum (see [Ml)- 

In this thesis, we have used methods from algebra, combinatorics, and polyhedral 
geometry to construct and enumerate magic squares and these methods are similar to 
the methods used by Stanley to enumerate semi-magic squares jHl]. Part of our work 
is based on generalizing his ideas. Polyhedral methods were also used by Halleck in 
his 2000 Ph.D thesis (see ^7\) and Beck et al. in [T2] to enumerate various kinds of 
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Figure 1.4: The pandiagonals. 

magic squares. Our techniques are different and are more algebraic in ffavor. 

The basic idea is to consider the defining hnear equations of magic squares. These 
equations, together with the nonnegativity requirement on the entries, imply that 
the set of magic squares becomes the set of integral points inside a pointed polyhedral 
cone (see or ^5)- The minimal Hilbert basis of this cone is defined to be the 
smallest finite set S of integral points with the property that any integral point 
can be expressed as a linear combination with nonnegative integer coefficients of the 
elements of S jHI]. For example, the minimal Hilbert basis of the 3x3 magic squares 
is given in Figure 11.51 and a Hilbert basis construction of the Loh-shu magic square 
is given in Figure 11.61 

102 201 021 120 111 

210 012 210 012 111 

021 120 102 201 111 

Figure 1.5: The minimal Hilbert Basis of 3 x 3 Magic squares. 



4 9 2 
3 5 7 
8 1 6 

Figure 1.6: A Hilbert basis construction of the Loh-Shu magic square. 

We map magic squares to monomials in a polynomial ring to enumerate them and 
derive the formula for the number of 4 x 4 magic squares of magic sum s (also derived 
simultaneously by Beck et al. \l2i using different techniques). 
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Theorem 1.1.1. The number 0/4x4 magic squares with magic sum s 
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480 



S 



,7 
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when s is even. 



We will describe this method in detail in the next section (also, see |2] and 0). 

1.2 Generating and enumerating lattice points inside poly- 
hedral cones. 

In this section we illustrate the algebraic techniques of generating and enumerating 
lattice points inside polyhedral cones by applying the method to construct and enu- 
merate magic squares. As an example, we walk through the details of the proof of 
Theorem ll.l.ll 

For these purposes we regard n x n magic squares as either n x n matrices or 

2 

vectors in and apply the normal algebraic operations to them. We also consider 
the entries of an n x n magic square as variables yij (1 < i, j < n). If we set the first 
row sum equal to all other mandatory sums, then magic squares become nonnegative 
integral solutions to a system of linear equations Ay = 0, where A is an {2n + 1) x 
matrix each of whose entries is 0, 1, or -1. 

For example, the equations defining 3x3 magic squares are: 
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yil + yi2 + 1/13 = Z/21 + 1/22 + Z/23 

yil + Z/12 + 1/13 = Z/31 + 1/32 + 1/33 

1/11 + 1/12 + 1/13 = 1/11 + 1/21 + 1/31 

1/11 + 1/12 + 1/13 = 1/12 + 1/22 + 1/32 

1/11 + 1/12 + 1/13 = 1/13 + 1/23 + 1/33 

1/11 + 1/12 + 1/13 = 1/11 + 1/22 + 1/33 

1/11 + 1/12 + 1/13 = 1/13 + 1/22 + 1/31 

Therefore, 3x3 magic squares are nonnegative integer solutions to the system of 
equations Ay = where: 



A 



1 1 1 

1 1 1 

1 1 

1 1 
1 1 

1 1 

1 1 



-1 


-1 







-1 




-1 



-1 
-1 



-1 





-1 








-1 -1 

-1 

-1 





-1 





-1 




-1 
-1 





and y 



1/11 

1/12 
1/13 
1/21 
1/22 
1/23 
1/31 
1/32 
1/33 



A nonempty set C of points in M" is a cone if au + bv belongs to C whenever u 
and V are elements of C and a and b are nonnegative real numbers. A cone is pointed 
if the origin is its only vertex (or minimal face; see [SI])- A cone C is polyhedral if 
C = {y : Ay < 0} for some matrix A, i.e, if C is the intersection of finitely many 
half-spaces. If, in addition, the entries of the matrix A are rational numbers, then 
C is called a rational polyhedral cone. A point y in the cone C is called an integral 
point if all its coordinates are integers. 

It is easy to verify that the sum of two magic squares is a magic square and 
that nonnegative integer multiples of magic squares are magic squares. Therefore, 
the set of magic squares is the set of all integral points inside a polyhedral cone 
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Cm„ = {y '■ Ay = 0,y > 0} in , where A is the coefficient matrix of the defining 
linear system of equations. Observe that Ca/„ is a pointed cone. 

In the example of 4 x 4 magic squares, there are three linear relations equating 
the first row sum to all other row sums and four more equating the first row sum 
to column sums. Similarly, equating the two diagonal sums to the first row sum 
generates two more linear equations. Thus, there are a total of 9 linear equations 
that define the cone of 4 x 4 magic squares. The coefficient matrix A has rank 8 and 
therefore the cone Cm4, of 4 x 4 magic squares has dimension 16 — 8 = 8 (see [SI])- 

In 1979, Giles and PuUeyblank introduced the notion of a Hilbert basis of a cone 
For a given cone C, its set Sc = C nZ^ of integral points is called the semigroup 
of the cone C. 

Definition 1.2.1. A Hilbert basis for a cone C is a finite set of points HB{C) in its 
semigroup Sc such that each element of Sc is a linear combination of elements from 
HB{C) with nonnegative integer coefficients. 

For example, the integral points inside and on the boundary of the parallelepiped 
in with vertices (0, 0), (3, 2), (1, 3) and (4, 5) in Figure IT77I form a Hilbert basis of 
the cone generated by the vectors (1, 3) and (3, 2). 




Figure 1.7: A Hilbert Basis of a two dimensional cone. 

We recall an important fact about Hilbert bases [SU Theorem 16.4]: 

Theorem 1.2.1. Each rational polyhedral cone C is generated by a Hilbert basis. 
If C is pointed, then there is a unique minimal integral Hilbert basis generating C 
(minimal relative to taking subsets). 
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We present a proof of Theorem 11.2.11 in Appendix ^ The minimal Hilbert basis 
of a pointed cone is unique and henceforth, when we say the Hilbert basis, we mean 
the minimal Hilbert basis. An integral point of a cone C is irreducible if it is not a 
linear combination with integer coefficients of other integral points. All the elements 
of the minimal Hilbert basis are irreducible jSHI • Since magic squares are integral 
points inside a cone. Theorem 11.2.11 implies that every magic square is a nonnegative 
integer linear combination of irreducible magic squares. 

The minimal Hilbert basis of the polyhedral cone of 4 x 4 magic squares is given in 
Figure HTHl and was computed using the software MLP (now called 4ti2) (see |SHI; soft- 
ware implementation 4ti2 is available from http:/ /www. 4ti2.de). In fact, all Hilbert 
basis calculations in this thesis are done with the software 4ti2. 

Hilbert basis constructions of the Jaina, and the Diirer magic squares are given in 
Figures 11.91 and II.IUI respectively. 
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0100 1000 1000 0101 0011 

0010 0010 0001 0020 1010 
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Figure 1.8: The minimal Hilbert Basis of 4 x 4 Magic squares. 



Different combinations of the elements of a Hilbert basis sometimes produce the 
same magic square. Figures 11.91 and 11.111 exhibit two different Hilbert basis con- 
structions of the Jaina magic square. This is due to algebraic dependencies among 
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Figure 1.9: A Hilbert basis construction of the Jaina magic square. 
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Figure 1.10: A Hilbert basis construction of Diirer's magic square. 



the elements of the Hilbert basis. Repetitions have to be avoided when count- 
ing squares, a problem that we solve by using a little bit of commutative algebra. 
Let HB{Cm„) = {hi, h2, . ■ ■ hr} be a Hilbert basis for the cone of n x n magic 
squares. Denote the entries of the square hp by yfj, and let k be any field. Let 
be the unique ring homomorphism between the polynomial rings k[xi,X2, ■ ■ ■ ,Xr] 

and k[tu,ti2, . . . ,tin,t21,t22, ■ ■ ■t2n, ■ ■ ■ ,tnl,tn2, ■ ■ ■ ,tnn] SUch that (p{Xp) = t^" , the 

monomial defined by 

i,j=l,...,n 

Monomials in k[xi,X2, . . . ,Xr] correspond to magic squares under this map, and 
multiplication of monomials corresponds to addition of magic squares. For example, 
the monomial x^xf^^ corresponds to the magic square 5hi + 200/i3. Different combi- 
nations of Hilbert basis elements that give rise to the same magic square can then 
be represented as polynomial equations. Thus, from the two different Hilbert basis 
constructions of the Jaina magic square represented in Figures 1131 and H . 1 11 we learn 
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Figure 1.11: Another Hilbert basis construction of the Jaina magic square. 

that 

/il + 4 ■ /i3 + 2 • M + 8 ■ /i5 + 3 • /i6 + 12 ■ /i7 + 4 ■ /i8 = 
h3 + 8-h5 + h6 + ll-h7 + h8 + hU + 2-hl5 + 2-hl7 + h20 
In k\ Xi, X2, ■ ■ ■ , Xr], this algebraic dependency of Hilbert basis elements translates to 

*^1*^3'^4*^5*^5'^7 "'^g iX'3*^5*^D'^7 *^o*^ 14*^ 15 11 -^0 — ^ • 

Consider the set of all polynomials in k[xi,X2, ■ ■ ■ , Xr] that are mapped to the zero 

polynomial under (p. This set, which corresponds to all the algebraic dependencies 

of Hilbert basis elements, forms an ideal in k[xi,X2, ■ ■ ■ ,Xr], an ideal known as the 

toric ideal of HB{Cm„) (see [E], or jSH] for details about toric ideals). If we 

denote the toric ideal as lHB(CM„)y then the monomials in the quotient ring Rcm„ = 

k[xi, X2, ■ ■ ■ , Xr]/ Ihb{Cm„) one-to-one correspondence with magic squares. 

For example, in the case of 3 x 3 magic squares, there are 5 Hilbert basis elements 
(see Figure IT75j) and hence there are 5 variables xi, X2, 0:3, X4, Xs which gets mapped 
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by as follows: 



Xi 



1 2 

2 1 
2 1 



X2 



2 1 
1 2 
12 



tll^t\it22t23^t3lt32^ 



X3 



2 1 
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tl2^tl3t21^t22t3lt33^ 



12 
1 2 
2 1 



t\ltl2^t22t23^t31^i'i'i 



X5 



1 1 
1 1 
1 1 



tllt\2tl3t2lt22t23t3lt32t33 



We use CoCoA to compute the toric ideal 



IhB{Cm^) = i^l^i - xhx2X3 - X1X4). 

Thus 

_ Q[xi,X2,X3,X4,X5] 

^ {XIX4 — X5, X2X3 — XIX4) 
Let Rcm„ (s) be the set of all homogeneous polynomials of degree s in the ring Rcm„ ■ 

Then Rc,^,,^{s) is a fc- vector space, and -Rc'^^^(O) = k. The dimension dimfc(i?c'A/„ (-s)) 

of -Rca/„ (■5) is precisely the number of monomials of degree s in -Rcm„ • 

Recall that if the variables Xi of a polynomial ring k[xi^X2i ■ ■ ■ iXr] are assigned 

nonnegative weights Wj, then the weighted degree of a monomial ■ ■ ■ x"'' is X]I=i ^i' 

Wi (see jS]). Therefore, if we take the weight of the variable Xi to be the magic sum of 

the corresponding Hilbert basis element /ij, then dimfc(i?(7j^^^ (s)) is exactly the number 
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of magic squares of magic sum s. For example, in the case of 3 x 3 magic squares, 
because all the elements of the Hilbert basis have sum 3, all the variables are assigned 
degree 3. 

Since Rcm„ a graded fc-algebra, it can be decomposed into a direct sum of its 
graded components Rcm„ = ® Rcm„{s) (see j2] or |H|). Consider a finitely generated 
graded /c-algebra i?CM„ = 0^Cm„(-5)- The function H{Rc^j^,s) = dim^ (/?c'm„ (-s)) is 
the Hilbert function of Rcm„ Hilbert- Poincare series of -Rcm„ is the formal 

power series 

oo 

We can now deduce the following lemma. 

Lemma 1.2.1. Let the weight of a variable Xi in the ring R = k[xi,X2, ...,Xr] be the 
magic sum of the corresponding element of the Hilbert basis hi. With this grading of 
degrees on the monomials of R, the number of distinct magic squares of magic sum s 
is given by the value of the Hilbert function H{Rcj,j^, s). 

We record here a version of the Hilbert-Serre Theorem. A proof of the Hilbert-Serre 
theorem in all its generality is presented in Appendix 1X1 

Theorem 1.2.2 (Theorem 11.1 [H]). Let k be a field and R := k[xi, X2, Xr] be 
a graded Noetherian ring, let xi,X2, ■■■,Xr be homogeneous of degrees > 0. Let M be 
a finitely generated R-module. Then the Hilbert Poincare series of M, HM(t) is a 
rational function of the form: 

where p{t) G 

By invoking the Hilbert-Serre theorem, we conclude that the Hilbert-Poincare 
series is a rational function of the form Hrc^^ {t) = p{t)/^i=i{^ ~ t'^'^^^'), where p{t) 
belongs to Z[t]. We refer the reader to jS], or for information about 
Hilbert-Poincare series. 

We can also arrive at the conclusion that H^^^^ (t) is a rational function by study- 
ing rational polytopes. A polytope V is called rational if each vertex of V has rational 
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coordinates. The dilation of a polytope V by an integer s is defined to be the polytope 
sV = {sa : a G V} (see Figure ITTT^ for an example). 



(2.4) (4.4) 



(1.2).,.. (2,2].. 



.■ (2,2) ■ 



(4,2) 



Figure 1.12: Dilation of a polytope. 



Let i(V, s) denote the number of integer points inside the polytope sV. If a G Q™", 
let den a be the least positive integer q such that qa G Z™. 

Theorem 1.2.3 (Theorem 4.6.25 p3j). Let V be a rational convex polytope of 
dimension d m M"* with vertex set V . Set F{V,t) = 1 + Xln>i^('^' '^)'^'^- ThenF{V,t) 
is a rational function, which can be written with denominator Y[aev(^ ~ fden ay 

To extract explicit formulas from the generating function we need to define the 
concept of quasi-polynomials. 

Definition 1.2.2. A function f : N ^ C is a quasi-polynomial if there exists an 
integer N > and polynomials fo, fi, fd such that 

f{n) = fi{n) if n = i{modN). 

The integer N is called a quasi-period of f . 

For example, the formula for the number of 3 x 3 magic squares of magic sum 
s is a quasi-polynomial with quasi-period 3 (see Section ITTT|) . We now state some 
properties of quasi-polynomials. 

Proposition 1.2.1 (Corollary 4.3.1 [53j). The following conditions on a function 
/ : N 1-^ C and integer N > are equivalent: 

1. f is a quasi-polynomial of quasi-period N . 
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,n 



P{x)_ 



Q(x) 



where P{x) and Q{x) G C[x], every zero a of Q{x) satisfies = 1 (provided 
P{x)/Q{x) has been reduced to lowest terms) and deg P < deg Q. 

3. For all n>0, 



where each Pi is a polynomial function of n and each 7^ satisfies 7^ = 1. The 
degree of Pi{n) is one less than the multiplicity of the root 7^^""^ in Q{x) provided 
P{x)/Q{x) has been reduced to lowest terms. 

Theorem ll.2.3l together with Proposition ll.2.l1 implv that i{V, s) is a quasi-polynomial 
and is generally called the Ehrhart quasi-polynomial ofV. A polytope is called an in- 
tegral polytope when all its vertices have integral coordinates. i{'P, s) is a polynomial 
if V is an integral polytope [SB]. 

The convex hull of n x n matrices, with nonnegative real entries, such that all the 
row sums, the column sums, and the diagonal sums equal 1, is called the polytope of 
stochastic magic squares. Then, clearly M„(s) is the Ehrhart quasi-polynomial of the 
polytope of stochastic magic squares. Therefore, by Theorem ll.2.3t we get, again, 
that Hj}„ it) is a rational function. 



Coming back to the case of 4 x 4 magic squares, we used the program CoCoA (see 
PT] : CoCoA software is available from http://cocoa.dima.unige.it) to compute the 
Hilbert-Poincare series YlT=o Mi{.s)t^ and obtained 



1 + 8t + 48^2 + 200t3 + 675t^ + 1904^^ + 4736^*^ + 10608t^ + 21925^^ + . . . 

Recall that the coefficient of is the number of magic squares of magic sum s. This 
information along with Proposition 11.2.11 enables us to recover the Hilbert function 
M4(s) in Theorem 1 1 . 1 . 1 1 from the Hilbert-Poincare series by interpolation. We provide 



' '1=1 



n 




s _ f^+4f"+18t'^+36fH50f''+36f^+18t^+4f+l 
(l-t)4{l-t2)4 
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basic algorithms to compute Hilbert bases, toric ideals, and Hilbert-Poincare series 
in Appendix iBl 

We can also study magic figures of higher dimensions with our methods. Similar 
enumerative results for magic cubes and further properties of the associated cones are 
discussed in Chapter |21 See Figure 11.131 for a Hilbert basis construction of a magic 
cube. 
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Figure 1.13: A Hilbert basis construction of a 3 x 3 x 3 magic cube. 



In the next section we apply these methods to construct and enumerate pandiag- 
onal magic squares. 

1.3 Pandiagonal magic squares. 

We continue the the study of pandiagonal magic squares started in [UISII. Here we 
investigate the Hilbert bases, and recompute the formulas of Halleck (computed using 
a different method, see |SZl Chapters 8 and 10]). 

The convex hull of n x n matrices, with nonnegative real entries, such that all the 
row sums, the column sums, and the pandiagonal sums equal 1, is called the polytope 
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of panstochastic magic squares. The integrality of the polytope of panstochastic magic 
squares was fully solved in pp. 

Let us denote by MP„(s) the number of n x n pandiagonal magic squares with 
magic sum s. Halleck computed the dimension of the cone to be {n — 2)^ for 
odd n and {n — 2)^ + 1 for even n (degree of the quasipolynomial MP„(s) is one less 
than these). For the 4x4 pandiagonal magic squares a fast calculation corroborates 
that there are 8 elements in the Hilbert basis (see Figure IT. 14j) . In his investigations, 
Halleck |37i identified a much larger generating set. Recall that the Jaina magic 
square is also a pandiagonal magic square. A pandiagonal Hilbert basis construction 
of the Jaina magic square is given in Figure IT. 151 



1 



1 1 

1 

1 1 



1 1 



1 
1 1 
1 



1 1 
1 1 



1 

1 1 



Figure 1.14: Hilbert basis of the 4x4 Pandiagonal magic squares. 



+8 



1 

1 



+2 



I I 



1 I 





+5 



Figure 1.15: A construction of the Jaina magic square with the Hilbert basis of pandiagonal magic 
squares. 



We verify that the 5x5 pandiagonal magic squares have indeed a polynomial 
counting formula. This case requires in fact no calculations thanks to earlier work by 
PP who proved that for = 5 the only pandiagonal rays are precisely the pandiagonal 
permutation matrices. It is easy to see that only 10 of the 120 permutation matrices 
of order 5 are pandiagonal. A 4ti2 computation shows that the set of pandiagonal 
permutation matrices is also the Hilbert basis (see Figure IT. l(i|) . 

We calculate the formulas stated in Theorem 11.3.11 using CoCoA: 
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Figure 1.16: Hilbcrt basis of the 5x5 Pandiagonal magic squares. 



Theorem 1.3.1. 



^(s^ + 4s + 12)(s + 2)2 if 2 divides s, 
otherwise. 



MP.(s) 



8064 



[s + 4)(s + 3)(s + 2)(s + l){s^ + 5s + 8){s^ + 5s + 42). 



In the next section, we apply our methods to study Frankhn squares which are 
more complex than the squares we have seen so far. 



1.4 Franklin Squares. 

The well-known squares Fl and F3, as well as the less familiar F2, that appear in 
Figure IT. 171 were constructed by Benjamin Franklin (see and |45]). In a letter to 
Peter Collinson (see |7j) he describes the properties of the 8x8 square Fl as follows: 

1. The entries of every row and column add to a common sum called the magic 
sum. 

2. In every half-row and half-column the entries add to half the magic sum. 

3. The entries of the main bent diagonals (see Figure IT. 19|) and all the bent diago- 
nals parallel to it (see Figure TTTIU^ add to the magic sum. 

4. The four corner entries together with the four middle entries add to the magic 
sum. 



CHAPTER 1. Introduction 



19 



52 


61 


4 


13 


20 


29 


36 


45 


14 


3 


62 


51 


46 


35 


30 


19 


53 


60 


5 


12 


21 


28 


37 


44 


11 


6 


59 


54 


43 


38 


27 


22 


55 


58 


7 


10 


23 


26 


39 


42 


9 


8 


57 


56 


41 


40 


25 


24 


50 


63 


2 


15 


18 


31 


34 


47 


16 


1 


64 


49 


48 


33 


32 


17 



200 217 232 249 8 



72 89 104 121 136 



Fl 



17 


47 


30 


36 


21 


43 


26 


40 


32 


34 


19 


45 


28 


38 


23 


41 


33 


31 


46 


20 


37 


27 


42 


24 


48 


18 


35 


29 


44 


22 


39 


25 


49 


15 


62 


4 


53 


11 


58 


8 


64 


2 


51 


13 


60 


6 


55 


9 


1 


63 


14 


52 


5 


59 


10 


56 


16 


50 


3 


61 


12 


54 


7 


57 



F2 



58 39 26 

i 219 230 

60 37 28 

201 216 233 

55 42 23 



7 250 

251 6 

5 252 

248 9 

10 247 



186 167 154 



70 91 102 



188 165 156 



105 
151 



135 122 

123 134 

133 124 

120 137 

138 119 



203 214 235 246 11 



53 44 21 



205 212 237 244 13 



75 86 107 118 139 



181 172 149 140 117 



77 84 109 116 141 



51 46 19 

207 210 239 

49 48 17 

196 221 228 

62 35 30 



14 243 

242 15 

16 241 

253 4 

3 254 



240 
29 
227 



209 
36 
222 



179 174 147 



79 82 II 



177 176 145 



93 
163 



194 223 226 255 2 



64 33 32 



142 115 

114 143 

144 113 

125 132 

131 126 

127 130 



164 

94 



192 161 160 129 128 



F3 



Figure 1.17: Squares constructed by Benjamin Franklin. 



Henceforth, when we say row sum, column sum, bent diagonal sum, and so forth, 
we mean that we are adding the entries in the corresponding configurations. Franklin 
mentions that the square Fl has five other curious properties but fails to list them. 
He also says, in the same letter, that the 16 x 16 square F3 has all the properties of 
the 8x8 square, but that in addition, every 4x4 subsquare adds to the common 
magic sum. More is true about this square F3. Observe that every 2x2 subsquare in 
F3 adds to one-fourth the magic sum. The 8x8 squares have magic sum 260 while 
the 16 X 16 square has magic sum 2056. For a detailed study of these three "Franklin" 
squares, see j7j, jlH], or 

We define 8x8 Franklin squares to be squares with nonnegative integer entries 
that have the properties (1) - (4) listed by Benjamin Franklin and the additional 
property that every 2x2 subsquare adds to one-half the magic sum (see Figure IT.18p . 
The 8x8 squares constructed by Franklin have this extra property (this might be 
one of the unstated curious properties to which Franklin was alluding in his letter). 
It is worth noticing that the fourth property listed by Benjamin Franklin becomes 
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: Magic sum 



: Magic sum 



: iialf the Magic sum 



: Magic sum 



: Magic sum = Magic sum 



: half the Magic sum = Magic sum 



: half the Magic sum 



: Magic sum 



Figure 1.18: Defining properties of the 8x8 Franklin squares 0- 



redundant with the assumption of this additional property. 

Similarly, we define 16 x 16 Franklin squares to be 16 x 16 squares that have non- 
negative integer entries with the property that all rows, columns, and bent diagonals 
add to the magic sum, the half-rows and half-columns add to one-half the magic sum, 
and the 2x2 subsquares add to one-fourth the magic sum. The 2x2 subsquare 
property implies that every 4x4 subsquare adds to the common magic sum. 

The property of the 2x2 subsquares adding to a common sum and the property 
of bent diagonals adding to the magic sum are "continuous properties." By this we 
mean that, if we imagine the square as the surface of a torus (i.e., if we glue opposite 
sides of the square together), then the bent diagonals and the 2x2 subsquares can 
be translated without effect on the corresponding sums (see Figure [1201) • 

When the entries of an n x n Franklin square (n = 8 or n = 16) are 1,2,3,..., n^, 
it is called a natural Franklin square. Observe that the squares in Figure 11.171 are 
natural Franklin squares. Very few such squares are known, for the simple reason 
that squares of this type are difficult to construct. Many authors have looked at the 
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Figure 1.19: The four main bent diagonals |45| . 






50: 63! 2 


15 


18 


31 


34: 47:50 


63 


2 


16 1 


64 


49 


48 


33 


32 


n: 16 


1 


64 


36 


45 


52 


61 


4 


13 


20 


29 


36 


45 


52 


61 


4 


30 


19 


14 


3 


62 


51 


46 


35 


30 


19 


14 


3 


62 


37 


44 


53 


60 


5 


12 


21 


28 


37 


44 


53 


60 


5 


27 


22 


11 


6 


59 


54 


43 


38 


27 


22 


11 


6 


59 


39 


42 


55 


58 


7 


10 


23 


26 


39 


42 


55 


58 


7 


25 


24 


9 


8 


57 


56 


41 


40 


25 


24 


9 


8 


57 


34 


47 


50 


63 


2 


15 


18 


31 


34 


47 


50 


63 


2 


32 
36 


17 
45 


16 


1 


64 


49 


48 


33 


32 


17 


16 


1 


64 


52 


61 


4 


13 20 


29 


36 


45 


52 


61 


4 




14: 3 


62 


51 : 46 


35 


30 


19 


14 


3 


62 



Figure 1.20: Continuous properties of Franklin squares. 



problem of constructing them (see, for example, (?], and jlTj and the references 
therein). 

Our method is a new way to construct Fl, F2, and F3. We are also able to 
construct new natural Franklin squares, not isomorphic to the ones previously known. 
A permutation of the entries of a Franklin square is a symmetry operation if it maps 
the set of all Franklin squares to itself, and two squares are called isomorphic if it is 
possible to transform one to the other by applying symmetry operations. We start 
by describing several symmetries of the Franklin squares. 

Theorem 1.4.1. The following operations on n x n Franklin squares, where n = 8 
or n = 16, are symmetry operations: rotating the square by 90 degrees; reflecting the 
square across one of its edges; transposing the square; interchanging alternate columns 
(respectively, rows) i and i + 2, where 1 < i < {n/2) — 2 or {n/2) + 1 < i < n — 2; 
interchanging the first n/2 columns (respectively, rows) of the square with the last n/2 
columns (respectively, rows) simultaneously; interchanging all the adjacent columns 
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(respectively, rows) i and i + 1 = 1, 3, 5, . . . n — 1) of the square simultaneously. 

The following operations are additional symmetry operations on 16 x 16 Franklin 
squares: interchanging columns (respectively, rows) i and i + A, where 1 < i < 4 or 
9 < z < 12. 

The proof of Theorem 11.4.11 is presented in Chapter El 
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To get this square 
from F2 interchange 
rows I and 3, 
rows 2 and 4, 
rows 5 and 7, 
rows 6 and 8, 
columns 1 and 3, 
columns 2 and 4, 
columns 5 and 7, 
and columns 6 and 8. 



To get this square 
from F3 interchange 
rows 1 and 3, 
rows 4 and 6, 
and rows 9 and 1 1 . 
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Figure 1.21: Constructing Franklin squares by row and column exchanges of Franklin squares. 



The following theorem addresses itself to the squares in Figures 11.171 and 11.221 

Theorem 1.4.2. The original Franklin squares Fl and F2 are not isomorphic. The 
squares Nl, N2, and N3 are nonisomorphic natural 8x8 Franklin squares that are 
not isomorphic to either Fl or F2. The square N4 is a natural 16 x 16 Franklin square 
that is not isomorphic to square F3. 

We give a proof of Theorem 11.4.21 in Chapter El We now enumerate Franklin 
squares with our methods. 

Theorem 1.4.3. Let Fs{s) denote the number o/8 x 8 Franklin squares with magic 
sum s. Then: 
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Figure 1.22: New natural Franklin squares constructed using Hilbert bases. 



Fais) = < 



_23 ,8 . 



359 „2 177967 



if s = 2 (mod 12) and s^2, 



23 9 , 23 „8 , 167 ,7 , 5 „6 , 581 „6 , 1823 ,4 , 6127 „3 

627056640 ^ 17418240 * ^ 6531840 ^ 15552 * ^ 186624 * ^ 77760 ^ 46360 



10741 2 _|_ 113443 , 3211 
20412 * 102060 2187 



23 9 , 

627056640 ^ 17418240 ' 6531840 



if s = i (mod 12), 



15552 " ' 933120 ' 



_5_ 2 _ 3967 „ _ 13 
378 * 10080 * 8 



if s = 6 ( mod 12 ), 



23 „9 , 23 „8 1 167 7 , 5 „6 , 581 „5 , 1823 „4 , 6127 „3 
627056640 17418240 6531840 15552 * 186624 77760 45360 



11189 J2 I 167203 



5771 
2187 



ifs = 8 (mod 12), 



23 „8 1 167 „7 I 5 „6 , 2419 „5 , 1013 „4 , 701 „3 
17418240 6531840 15552 * 933120 77760 22680 



10206 816480 17496 



if ss 10 (mod 12), 
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The 8x8 pandiagonal Franklin squares are the 8x8 Frankhn squares that have all 
the pandiagonals adding to the common magic sum (see Figure EU- Our techniques 
enable us to construct and count 8x8 pandiagonal Franklin squares as well: 

Theorem 1.4.4. Let PFs{s) denote the number of 8x8 pandiagonal Franklin squares 
with magic sum s. Then: 



PF.(s) = I 



2293760 ^ 71680 ^ 3840 ^ 320 ^40 ^15 ^ 420 ^ 105 ^ 

if s = ( mod 4 ), 



otherwise. 

We introduce magic graphs and explore its connections to symmetric magic squares 
in the next section. 



1.5 Magic labelings of graphs. 

Let G be a finite graph. A labeling of G is an assignment of a nonnegative integer 
to each edge of G. A magic labeling of magic sum r of G is a labeling such that 
for each vertex f of G the sum of the labels of all edges incident to v is the magic 
sum r (loops are counted as incident only once). Graphs with a magic labeling are 
also called magic graphs fsee 11.231 for an example of a magic labeling of the complete 
graph Kq). Magic graphs are also studied in great detail by Stanley and Stewart in 

ISSl, ISni, ISIl, and PI. 

We define a magic labeling of a digraph D of magic sum r to be an assignment of 
a nonnegative integer to each edge of D, such that for each vertex Vi of D, the sum 
of the labels of all edges with Vi as the initial vertex is r, and the sum of the labels of 
all edges with Vi as the terminal vertex is also r. Thus magic labelings of a digraph 
is a network flow, where the flow into and out of every vertex, is the magic sum of 
the labeling (see for details about flows). Interesting examples of magic digraphs 
are Cayley digraphs of flnite groups. Let G be a flnite group {gi,g2, ■ ■ ■ ,gn = I}- 
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1 




Figure 1.23: A magic labeling of the complete graph Kq of magic sum 40 |58| . 

The Cayley group digraph of G is a graphical representation of G: every element Qi 
of the group G corresponds to a vertex Vi {i = 1,2, ... ,n) and every pair of distinct 
vertices Vi,Vj is joined by an edge labeled with a where ga = djdi^ (see [40^ or [^T]). 
For example, the Cayley digraph for the permutation group 

^3 = {gi = (123), (?2 = (132), = (23), = (12), (75 = (13), (76 = /} 
is given in Figure I1.241 

Proposition 1.5.1. The Cayley digraph of a group of order n is a magic digraph 
with magic sum "■^'^~^'> . 

Proof. Let Cij denote an edge between the vertex Vi and the vertex Vj of the Cayley 
digraph such that Vi is the initial vertex and Vj is the terminal vertex. Let vi be a 
vertex of the Cayley digraph, and let a be an integer in the set {1,2, . . . ,n — 1}. 
Let gp = g^gi and let gq = giga- Then, the edges cip and Cgi are labeled by a. Also, 
gjg^^ = (7„ = / if and only ii i = j. Hence, a Cayley group digraph is a magic digraph 
with magic sum 1 + 2 + - — |-(n — 1) = (see also Chapter 8, Section 5 in |1T|). 

□ 

A digraph is called Eulerian if for each vertex v the indegree and the outdegree 
of V is the same. Therefore, Eulerian digraphs can also be studied as magic digraphs 
where all the edges are labeled by 1 (see jS] for the applications of Eulerian digraphs 
to digraph colorings). 
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Figure 1.24: Cayley digraph of the group 5*3 |41| . 

If we consider the labels of the edges of a magic graph G as variables, then the 
defining magic sum conditions are linear equations. Thus, like before, the set of magic 
labelings of G becomes the set of integral points inside a pointed polyhedral cone Cg 
(see Section ITT^ . Therefore, with our methods we can construct and enumerate magic 
labelings of graphs. For example, a Hilbert basis construction of a magic labeling of 
the complete graph Kq is given in Figure IT. 251 

Let Hg{t) denote the number of magic labelings of G of magic sum r. The generat- 
ing functions of Heir) in this thesis were computed using the software LattE (see P7j: 
software implementation LattE is available from http: / /www. math. ucdavis.edu/ latte) 
LattE was able to handle computations that CoCoA was not able to perform. LattE 
uses Barvinok's algorithm to compute generating functions which is different from 
the methods used by CoCoA (see |2Z|)- 

For example, let r„ denote the complete general graph on n vertices (i.e, the com- 
plete graph with one loop at every vertex). The formulas ifr„(^) for ^ = 3 and n = 4 
were computed by Carlitz [22 , and Hr^{r) was computed by Stanley ^56j. We use 
LattE to derive H-pg (r) : 
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Figure 1.25: A Hilbert basis construction of a magic labeling of the complete graph K^,. 



243653 15 i 243653 14 , 91173671 13 i 5954623 12 



1992646656000 ' 44281036800 ' 797058662400 ' 4087480320 



I 3895930519 11 , 21348281 10 i 1063362673 9 , 7132193 8 , 479710409 7 



306561024000 ' 265420800 ' 2786918400 ' 5160960 ' 124416000 



I 963567863 6 i 26240714351 5 i 39000163 4 , 1514268697 3 



116121600 ' 1916006400 ' 2280960 ' 96096000 



74169463^2 , 176711 



2 I VlKiill I -1 



7207200 ' 40040 



if 2 divides r, 



^re(r) = <^ 



243653 15 , 243653 14 , 91173671 13 , 5954623 12 



1992646656000' ' 44281036800' ' 797058662400' ' 4087480320' 



I 3895930519 ^11 _|_ 21348281 ^10 _|_ 1063362673 ^9 _|_ 7132193 ^8 _|_ 479710409 ^7 



306561024000 ' 265420800 ' 2786918400 ' 5160960 ' 124416000 



I 963567863 6 i 839695842607 5 i 9983039353 4 , 774706849739 3 
116121600 61312204800 583925760 49201152000 
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An n-matching of G is a magic labeling of G with magic sum at most n and the 
labels are from the set {0, 1, ... , n} (see jl2], chapter 6). A perfect matching of G is 
a 1-matching of G with magic sum 1. 

Proposition 1.5.2. The perfect matchings of G are the minimal Hilbert basis ele- 
ments of Gg of magic sum 1 and the number of perfect matchings of G is Hg{1)- 

Proof. Magic labelings of magic sum 1 always belong to the minimal Hilbert basis 
because they are irreducible. Therefore, perfect matchings belong to the minimal 
Hilbert basis because they have magic sum 1. Conversely, every magic labeling of 
magic sum 1 is a perfect matching. So we conclude that the perfect matchings of 
G are the minimal Hilbert basis elements of Gg of magic sum 1. The fact that the 
number of perfect matchings of G is Hg{^) follows by the definition of Hg{1). □ 

The perfect matchings of G can be found by computing a truncated Hilbert basis 
of magic sum 1 using 4ti2 (see ^S])- 

Hilbert bases can also be used to study factorizations of labeled graphs. We define 
Factors of a graph G with a labeling L to be labelings Li,i = 1, ... ,r of G such that 
-^(^) = Si=i Li{G), and if Li{ek) ^ for some edge Ck of G, then Lj{ek) = for all 
j ^ i. A decomposition of L into factors is called a factorization of G. An example 
of a graph factorization is given in Figure IT. 261 See Chapters 11 and 12 of |3T] for a 
detailed study of graph factorizations. 

Edye Label 




Figure 1.26: Graph Factorization. 

With our methods we can also construct and enumerate magic labelings of di- 
graphs. Let Ho{r) denote the number of magic labelings of a digraph D of magic 
sum r. We now connect the magic labelings of digraphs to magic labelings of bipartite 
graphs. 

Lemma 1.5.1. For every digraph D, there is a bipartite graph Go such that the 
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magic lahelings of D are in one-to-one correspondence with the magic lahelings of 
Gd- Moreover, the magic sums of the corresponding magic lahelings of D and Gn are 
also the same. 

Proof. Denote a directed edge of a digraph D with Vi as the initial vertex, and 
Vj as the terminal vertex, by Cij. Let L be a magic labeling of D of magic sum r. 
Consider a bipartite graph G/j in 2n vertices, where the vertices are partitioned into 
two sets A = {ai, . . . , a„} and B = {61, ... , such that there is an edge between 
Oi and bj, if and only if, there is an edge Cij in D. Consider a labeling Lq^ oi Gd 
such that the edge between the vertices and bj is labeled with L{eij). Observe that 
the sum of the labels of the edges incident to a, is the same as the sum of the labels 
of incoming edges at the vertex Vi of D. Also, the sum of the labels of edges at a 
vertex bj is the sum of the labels of outgoing edges at the vertex Vj of D. Since L 
is a magic labeling, it follows that Lg^ is a magic labeling of Gd with magic sum 
r. Going back-wards, consider a magic labeling L' of Gd- We label every edge Cij of 
D with the label of the edge between and bj of G/j to get a magic labeling of 
D. Observe that V and Ld have the same magic sum. Hence, there is a one-to-one 
correspondence between the magic labelings of D and the magic labelings of Gd. □ 

For example, the magic labelings of the Octahedral digraph with the given orien- 
tation Do in Figure IT. 271 are in one-to-one correspondence with the magic labelings 
of the bipartite graph Gdq- 




3 

Figure 1.27: The octahedral digraph with a given orientation Dq and its corresponding bipartite 
graph Gdo ■ 
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Similarly, for a bipartite graph B we can get a digraph such that the magic 
labelings of B are in one-to-one correspondence with the magic labelings oi Bq: let 
B be such that the vertices are partitioned into sets A = {ai,...,a„} and B = 
{bi, . . . ,bm}- Without loss of generality assume n > m. Then Bd is the digraph 
with n vertices such that there is an edge Cij in Ba ii and only if there is an edge 
between the vertices and hj in B. This correspondence enables us to generate and 
enumerate perfect matchings of bipartite graphs. 

Proposition 1.5.3. There is a one-to-one correspondence between the perfect match- 
ings of a bipartite graph B and the elements of the Hilbert basis ofCso- The number 
of perfect matchings of B is Hbjj{1)- 

We present a proof of Proposition ll.5.3l in ChapterHJ A graph G is called a positive 
graph if for any edge e of G there is a magic labeling L of G for which L(e) > j55j. 
Since edges of G that are always labeled zero for any magic labeling of G may be 
ignored to study magic labelings, we will concentrate on positive graphs in general. 
We use the following results by Stanley from |55| and [Snj to prove Theorems 11.5.41 
and 11.5.51 and Corollarv 11.5.3.11 

Theorem 1.5.1 (Theorem 1.1, [56j). Let G be a finite positive graph. Then either 
Hair) is the Kronecker delta 5or or else there exist polynomials /g(^) one? Jci^) such 
that Heir) = lair) + {-If J air) for all r G N. 

Theorem 1.5.2 (Theorem 1.2, [56j). Let G be a finite positive graph with at least 
one edge. The degree of Heir) is q — n + b, where q is the number of edges of G, n is 
the number of vertices, and b is the number of connected components of G which are 
bipartite. 

Theorem 1.5.3 (Theorem 1.2, | j55j ^ . Let G be a finite positive bipartite graph with 
at least one edge, then Heir) is a polynomial. 

We now conclude that H£,{r) is a polynomial for every digraph D. 
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Corollary 1.5.3.1. If D is a digraph, then Hij{r) is a polynomial of degree q — 2n + b, 
where q is the number of edges of D, n is the number of vertices, and b is the number 
of connected components of the bipartite graph Go. 

Proof. The one-to one correspondence between the magic labelings of D and the 
magic labelings of Gd, implies by Theorem 11.5.31 that H£){r) is a polynomial, and 
by Theorem 11.5.21 that the degree of Hu{r) is q — 2n + b, where b is the number of 
connected components of Go that are bipartite. □ 

Consider the polytope V := < b}. Let c be a nonzero vector, and let S = 

max {cxl^x < b}. The affine hyperplane {x\cx = 6} is called a supporting hyperplane 
of V. A subset F of P is called a face of P if F = P or if F is the intersection of V 
with a supporting hyperplane of V. Alternatively, F is a face of V if and only if F is 
nonempty and 

F = {xe V\A'x = b'} 

for some subsystem A'x < b' of Ax < b. See [HI] for basic definitions with regards to 
polytopes. 

Let fi, f2, . . . , fn denote the vertices of a graph G and let e^j, e^j, . . . , Cj^^ denote 
the edges of G that are incident to the vertex Vi of G. Consider the polytope 

mi 

VG = {LeCG^ M", ^ L(eiJ = 1; 2 = 1, . . . , n}. 

j=i 

We will refer to Vg as the polytope of magic labelings of G. Then, Hg{t) is the 
Ehrhart quasi-polynomial of Vg (see Section (T!^ . A face of Vg is a polytope of the 
form 

where Eq = {cjj, . . . , Cj^} is a subset of the set of edges of G. 

Theorem 1.5.4. Let G be a finite positive graph with at least one edge. Then the 
polytope of magic labelings of G, Vg is a rational polytope with dimension q — n + b, 
where q is the number of edges ofG, n is the number of vertices, and b is the number of 
connected components of G that are bipartite. The d- dimensional faces ofVG are the 
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d- dimensional polytopes of magic labelings of positive subgraphs of G with n vertices 
and at most n — b + d edges. 

We prove Theorem 11.5.41 in Chapter EJ Observe from Theorem 1 1 . 5 . 41 1 hat there is 
an edge between two vertices Vi and vj of Vg if and only if there is a graph with at 
most n — b + 1 edges, with magic labehngs Vi and Vj. The edge graph of Vr-^ is given 
in Figure IT. 281 Similarly, we can draw the face poset of Vg (see Figure IT. 291 for the 
face poset of Vr^)- 




Figure 1.28: The edge graph of Vr^- 



An n X n semi-magic square of magic sum r is an n x tt, matrix with nonnegative 
integer entries such that the entries of every row and column add to r. Doubly 
stochastic matrices are n x n matrices in such that their rows and columns add 
to 1. The set of all n x n doubly stochastic matrices form a polytope Bn, called the 
Birkhoff polytope. See ^U], ^H]; or [HI] for a detailed study of the Birkhoff polytope. 

A symmetric magic square is a semi-magic square that is also a symmetric matrix. 
Let Hn{r) denote the number of symmetric magic squares of magic sum r (see P^ . 

and jSni for the enumeration of symmetric magic squares). We define the poly- 
tope Sn ofnxn symmetric magic squares to be the convex hull of all real nonnegative 
n X n symmetric matrices such that the entries of each row (and therefore column) 
add to one. 

A one-to-one correspondence between symmetric magic squares M = [rriij] of magic 
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i\ m m m if 




Figure 1.29: The face poset oiVra 



sum r, and magic labelings of the graph of the same magic sum r was estabhshed 
in [HE]: let denote an edge between the vertex Vi and the vertex vj of r„. Label the 
edge Cij of r„ with rriij, then this labeling is a magic labeling of r„ with magic sum 
r. See Figure [00 for an example. Therefore, we get Vr„ = Sn and Hr„{r) = Hn{r). 



1 1 

7 5 




Figure 1.30: A magic labeling of Fs and its corresponding symmetric magic square. 



Corollary 1.5.4.1. The polytope of magic labelings of the complete general graph 
Vr„ is an n{n — l)/2 dimensional rational polytope with the following description 

Pr„ = {L = {L{e,,)eR^; L{eij) > 0;1 < t, j < n,t < j, 

Ej=i + EJ=i+i = 1 fori = l,...,n}. 

The d- dimensional faces of Vr„ are d- dimensional polytopes of magic labelings of 
positive graphs with n vertices and at most n + d edges. There are ('^"~^") faces of 
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Vr2n that are copies of the Birkhoff polytope B^. 

See Chapter 0] for the proof of Corollary I1.5.4.11 Again, as in the case of graphs, 
we define a polytope Vd of magic labelings of D: Let e^^, e^j, . . . , Cj^. denote the edges 
of D that have the vertex Vi as the initial vertex and let fi-^, fi^, . . . , fi^_ denote the 
edges of D for which the vertex Vi is the terminal vertex, then: 

rrii Si 

VD = {LeCD^ R', L{e^,) = J2 ^(4) = 1; ^ = 1, ■ ■ • , ri}. 

We define a digraph D to be a positive digraph if the corresponding bipartite graph 
is positive. 

Theorem 1.5.5. Let D be a positive digraph with at least one edge. Then, Vd is an 
integral polytope with dimension q — 2n + b, where q is the number of edges of D, n is 
the number of vertices, and b is the number of connected components of Gd that are 
bipartite. The d-dimensional faces of Vd o^^e the d-dimensional polytopes of magic 
labelings of positive subdigraphs of D with n vertices and at most 2n — b + d edges. 

See Chapter |3] for the proof of Theorem II. 5. 51 

Let n„ denote the complete digraph with n vertices, i.e, there is an edge from 
each vertex to every other, including the vertex itself (thereby creating a loop at 
every vertex), then Gn„ is the the complete bipartite graph Kn,n- We get a one- 
to-one correspondence between semi-magic squares M = [rriij] of magic sum r and 
magic labelings of n„ of the same magic sum r by labeling the edges Cij of n„ with 
rriij. This also implies that there is a one-to-one correspondence between semi-magic 
squares and magic labelings of Kn,n (this correspondence is also mentioned in 
and [SZ|)- See Figure IT3T] for an example. 

A good description of the faces of Birkhoff polytope is not known jUj. We can 
now give an explicit description of the faces of the Birkhoff polytope. 

Theorem 1.5.6. Vjin the Birkhoff polytope Bn. The d-dimensional faces of Bn 
are d-dimensional polytopes of magic labelings of positive digraphs with n vertices and 
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Figure 1.31: Two different graph labelings associated to a semi- magic square. 



at most 2n + d — 1 edges. The vertices ofVo, where D is a positive digraph, are 
permutation matrices. 

The proof of Theorem ll.5.6l is presented in Chapter^ See Figure IT.32l for the edge 
graph of i?3. Two faces of a polytope of magic labehngs of a graph (or a digraph) are 
said to be isomorphic faces if the subgraphs (subdigraphs, respectively) defining the 
faces are isomorphic. A set of faces is said to be a generating set of d-dimensional 
faces if every c?- dimensional face is isomorphic to one of the faces in the set. See 
Figures 11.331 11.341 11.351 and 11.361 for the generators of the edges, the two dimensional 
faces, the facets, and the Birkhoff polytope -B3, respectively (the numbers in the 
square brackets indicate the number of faces in the isomorphism class of the given 
face). 
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Figure 1.33: The generators of the edges of the Birkhoff Polytope B3. 



[3] [6] [6] [6] [6] [31 




Figure 1.34: The generators of the 2-dimensional faces of the Birkhoff Polytope S3. 



[6] 



[3] 





Figure 1.35: The generators of the facets of the Birkhoff Polytope B3. 
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Figure 1.36: The Birkhoff Polytope B3. 
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Chapter 2 



Magic Cubes 



They flash upon that inward eye 

Which is the bliss of solitude; 

And then my heart with pleasure flUs, 

And dances with the daffodils. 

- William Wordsworth. 
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Figure 2.1: A magic cube. 



A semi-magic hypercube is a (i-dimensional nxnx ■ ■ ■ xn array of n'^ non-negative 
integers, which sum up to the same number s for any hne parallel to some axis. A 
magic hypercube is a semi-magic cube that has the additional property that the sums 
of all the main diagonals, the 2'^^^ copies of the diagonal X2,2,...,25 • • • ,Xn,n,...,n 

under the symmetries of the d-cuhe, are also equal to the magic sum. For example, 
in a 2 X 2 X 2 cube there are 4 diagonals with sums + ^2,2,2 = 2:2,1,1 + 2:1,2,2 = 

2^1,1,2 + 2^2,2,1 = 2^1,2,1 + 2;2,1,2- 
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An example of a 3 x 3 x 3 magic cube is given in Figure 12. 1[ If we consider the 
entries of a magic cube to be variables, the defining equations form a linear system 
of equations and thus magic cubes are integral points inside a pointed polyhedral 
cone. Therefore, we can use the methods described in Section 11.21 to construct and 
enumerate magic cubes. Let G denote the group of rotations of a cube Two cubes 
are called isomorphic if we can get one from the other by using a series of rotations. 
A set of magic cubes are called generators of the Hilbert basis if every element of 
the Hilbert basis is isomorphic to one of the cubes in the set. The generators of the 
Hilbert basis of 3 x 3 x 3 magic cubes are given in Figure YI7^ fthe numbers in square 
brackets indicate the number of elements in the orbit of a generator under the action 
of G). There are 19 elements in the Hilbert basis and all of them have magic sum 
value of 3. An example of a Hilbert basis construction of a magic cube is given in 
Figure 11.131 
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Figure 2.2: The generators of the Hilbert basis of the 3x3x3 magic cube. 



Let MC„(s) denote the number of x n x magic cubes of magic sum s. We 
use the algorithm presented in Section [T!^ to compute the generating function for the 
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number of 3 x 3 x 3 magic cubes: 

= 1 + 19 + 121 16 + 439 19 + 1171 + 2581 1^^ + 4999 1^^ + 
and we derive 

Theorem 2.0.7. The number 0/3x3x3 magic cubes 



otherwise. 



Theorem 2.0.8. The number of n x n x n magic cubes of magic sum s, MCn{s) is 
a quasipolynomial of degree {n — 1)^ — 4 for n > 3,n ^ 4. For n = 4 it has degree 
(4-1)3-3 = 24. 

Proof. The function that counts magic cubes is a quasipolynomial whose degree is 
the same as the dimension of the cone of magic cubes minus one. For small values 
(e.g = 3, 4) we can directly compute this. We present an argument for its value for 
n > 4. Let B be the (3n^ + 4) x matrix with 0, 1 entries determining axial and 
diagonal sums. In this way we see that nxnxn magic cubes of magic sum s are the 
integer solutions of Bx = (s, s, . . . , s)"^ , x > 0. 

It is known that for semi-magic cubes the dimension is {n — 1)^ ^2], which means 
that the rank of the submatrix B' of B without the 4 rows that state diagonal sums 
is — (n — 1)3. It remains to be shown that the addition of the 4 sum constraints 
on the main diagonals to the defining equations of the nxnxn semi-magic cube 
increases the rank of the defining matrix B by exactly 4. 

Let us denote the n^ entries of the cube by Xi^i^i, . . . , Xn,n,n and consider the {n — 
1) X [n — 1) X [n — 1) sub-cube with entries Xi^i^i, . . . , Xn-i,n-i.n~i- For a semi-magic 
cube we have complete freedom to choose these (n — l)^ entries. The remaining entries 
of the nxnxn magic cube become known via the semi-magic cube equations, and 
all entries together form a semi-magic cube. For example: 

i=l •^l,n,n 2_/i=l 2_/j=l •^i,3,'^i •^n,n,n 2_/j=l Z_/j=l 2_/fc=l 
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However, for the magic cube, 4 more conditions have to be satisfied along the 
main diagonals. Employing the above semi-magic cube equations, we can rewrite 
these 4 equations for the main diagonals such that they involve only the variables 
Xi^i^i, . . . ,Xn-i,n-i,n-i- Thus, as we will see, the complete freedom of choosing val- 
ues for the variables xi^i^i, . . . ,Xn-i,n~i,n-i is restricted by 4 independent equations. 
Therefore the dimension of the kernel of B is reduced by 4. 

Let us consider the 3 equations in Xi^i^i, . . . , Xn-i,n-i,n-i corresponding to the main 
diagonals . . . , Xn,n,i, Xi^ri,i, Xn,i,n, and Xn,i,i, xi^n,n- They are linearly 

independent, since the variables Xn-i,n-i,i, Xn-i,i,n-i^ and xi^n-i,n-i appear in exactly 
one of these equations. The equation corresponding to the diagonal xi^i^i, . . . , Xn,n,n 
is linearly independent from the other 3, because, when rewritten in terms of only 
variables of the form Xij^k with 1 < i, j, k < n, it contains the variable ^2,2,3; which 
for n > 4 does not lie on a main diagonal and is therefore not involved in one of the 
other 3 equations. Therefore, for n > A the kernel of the matrix B has dimension 
(n — 1)^ — 4. This completes the proof. □ 

Similarly, we can construct and enumerate semi-magic cubes. Bona IT] had al- 
ready observed that a Hilbert basis must contain only elements of magic constant one 
and two. Here, we compute the 12 Hilbert basis elements of magic sum 1 and the 54 
elements of magic sum 2 using 4ti2. The generators of the Hilbert basis of 3 x 3 x 3 
semi-magic cubes are given in Figure 12.31 (the definition of a generating set of the 
Hilbert basis of semi-magic cubes is similar to the corresponding definition for magic 
cubes). 

Denote by SH^{s) the number of semi-magic d-dimensional hypercubes with n"^ 
entries. We use CoCoA to compute the generating function SH^{s) : 

V^oo qtt3(„\+s _ t^+5t'^+67t''+130t"'+242f*+130t^+67t^+5t+l 

= 1 + 12t + 132^2 + 847t3 + 3921t^ + U286t^ + 43687^^ + 116757t^ + . . . . 

In ^Tj, Bona presented a proof that the counting function of 3 x 3 x 3 semi- 
magic cubes is a quasi-polynomial of non-trivial period. We improve on his result by 
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Figure 2.3: The generators of the minimal Hilbert basis of the 3x3x3 semi-magic cube. 



computing an explicit formula. 

Theorem 2.0.9. The number o/3 x 3 x 3 semi-magic cubes of magic sum s, 

9 o8 I ^ „7 I ^ 6 I 297 5 I 1341 ,4 , 513 3 , 3653 2 i 627 , i 
2240 ~'~ 560 ~'~ 320 ~'~ 320 ~'~ 640 ~'~ 160 ~'~ 1120 ~'~ 280 ~'~ 



SHiis) 



if2\s, 



9 ,8 I ^ 7 I ^ 6 I 297 5 I 1341 ^4 , 513 3 , 3653 2 , 4071 , 47 
2240 ~'~ 560 ~'~ 320 ~'~ 320 ~'~ 640 ~'~ 160 ~'~ 1120 ~'~ 2240 ~'~ 128 

otherwise. 



The convex hull of all real nonnegative semi-magic cubes (of given size) all whose 
mandated sums equal 1 is called the polytope of stochastic semi-magic cubes. The 
polytope of 3 X 3 X 3 stochastic semi-magic cubes is actually not equal to the convex 
hull of integral semi-magic cubes. This is because the 54 elements of degree two in the 
Hilbert basis, when appropriately normalized, give rational stochastic matrices that 
are all vertices. In other words, the Birkhoff-von Neumann theorem jMJ page 108] 
about stochastic semi-magic matrices is false for 3x3x3 stochastic semi-magic cubes. 
We prove the following result about the number of vertices of stochastic semi-magic 
cubes. 
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Theorem 2.0.10. The number of vertices of the polytope of n x n x n stochastic 
semi-magic cubes is bounded below by (?7,!)^"/n"^. 

Proof. We exhibit a bijection between integral stochastic semi-magic cubes and nxn 
latin squares: Each 2-dimensional layer or slice of the integral stochastic cubes are 
permutation matrices (by Birkhoff-Von Neumann theorem), the different slices or 
layers cannot have overlapping entries else that would violate the fact that along a 
line the sum of the entries equals one. Thus make the permutation coming from the 
first slice be the first row of the latin square, the second slice permutation gives the 
second row of the latin square, etc. From well-known bounds for latin squares we 
obtain the lower bound (see Theorem 17.2 in [HS])- D 
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Chapter 3 

Franklin Squares 

If back we look on ancient Sages Schemes, 
They seem ridiculous as Childrens Dreams 
- Benjamin Franklin. 

3.1 All about 8x8 Franklin squares. 

Like in the case of magic squares, we consider the entries of an n x n Franklin square 
as variables yij {1 < i, j < n) and set the first row sum equal to all other mandatory 
sums. Thus, Franklin squares become nonnegative integral solutions to a system of 
linear equations Ay — 0, where A is an (n^ + 8n — 1) x matrix each of whose entries 
is 0, 1, or -1. 

In the case of the 8x8 Franklin squares, there are seven linear relations equating 
the first row sum to all other row sums and eight more equating the first row sum to 
column sums. Similarly, equating the eight half-row sums and the eight half-column 
sums to the first row sum generates sixteen linear equations. Equating the four 
sets of parallel bent diagonal sums to the first row sum produces another thirty-two 
equations. We obtain a further sixty-four equations by setting all the 2x2 subsquare 
sums equal to the first row sum. Thus, there are a total of 127 linear equations that 
define the cone of 8 x 8 Franklin squares. The coefficient matrix A has rank 54 and 
therefore the cone C of 8 x 8 Franklin squares has dimension 10. 
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Let {ri,rj) denote the operator that acts on the space of n x n matrices by in- 
terchanging rows i and j of each matrix, and let {ci,Cj) signify the corresponding 
operator on columns. 

Consider the group G of symmetry operations of 8 x 8 Franklin squares (see Lemma 
I3.3.1|) : G is generated by the set 

{(ci, cs), (C5, Ct), (c2, C4), (C6, Cg), (ri, Tg), (rg, tt), (rs, r4), (re, rg)}. 

The Hilbert basis of the polyhedral cone of 8 x 8 Franklin squares is generated by 
the action of the group G on the three squares Tl, T2, and T3 in Figure EIH and 
their counterclockwise rotations through 90 degree angles. Not all squares generated 
by these operations are distinct. Let R denote the operation of rotating a square 90 
degrees in the counterclockwise direction. Observe that R^-Tl is the same as Tl and 
R^-Tl coincides with R-Tl. Similarly, R^-T2 is just T2, and R^-T2 is the same as 
R-T2. Also Tl and i?-Tl are invariant under the action of the group G. Therefore 
the Hilbert basis of the polyhedral cone of 8 x 8 Franklin squares consists of the 
ninety-eight Franklin squares: Tl and -R-Tl; the thirty-two squares generated by the 
action of G on T2 and -R-T2; the sixty-four squares generated by the action of G on 
T3 and its three rotations -R-T3, i?^-T3, and R^-T3. 
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Figure 3.1: Generators of the Hilbert basis of 8 x 8 Franklin squares. 



The Hilbert basis constructions of the Franklin squares F2, Nl, N2, Fl, and N3 
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F2 = 5 ■ hi + 16 ■ h2 + 4h3 + 3 ■ h4 + 2 ■ h5 + h6 + +32 ■ h7 + 2 ■ h8; 
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Nl = 5 ■ hi + 2 ■ h2 + 4 ■ h3 + 3 ■ h4 + 2 ■ h5 + h6 + 32 ■ h7 + 16 ■ h8; 

N2 = 5-hl + 16-h2 + 4-(c5,C7)-(c6,C8)-h3 + 3-h4 + 2-h5 + h6 + 32-h7 
+2 ■ h8; 

Fl = 2-hl + 14-(cl,c3)-hl + h2 + (r6,r8)-h2 + 3-(rl,r3)-(r6,r8)-h2 

+30 ■ h6 + 2 ■ (c5, c7) ■ h6 + 6 ■ h8 + 3 ■ (c5, c7) ■ T3 + (r2, r4) ■ (c5, c7) ■ T3; 

N3 = 2 ■ hi + 6 ■ (cl,c3) ■ hi + h2 + (r6,r8) ■h2 + 3 ■ (rl,r3) • (r6,r8) ■h2 

+30 ■ h6 + 2 ■ (c5, c7) ■ h6 + 14 ■ h8 + 3 ■ (c5, c7) ■ T3 + (r2, r4) • (c5, c7) ■ T3. 
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Figure 3.2: Constructing Benjamin Franklin's 8x8 square F2. 
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Similarly, let 



1 — 
J- — 


yr3, r^} 




V 11? ' 13j V 12, ' 14j "J^, 8^ — 


01, go — rt ■ gi. 


zL — 
^ — 


i <Y> ^ ly ^ \ 




V 9, ' I3J V 127 ' 16j 'J'^j go — 


LicLllbpUoc Ui OZr, 


6 = 


(cg, Cis) 


■ (cio,Ci4 


) ■ (C11,C15) ■ (ci2,Ci6) ■ g5. 




7 = 






■ (^^12,^6) ■ S2, 




8 = 


S3, g9 


= (n,r5) 


■ {rz.rj) ■ (rio,ri4) ■ (ri2,ri6) ■ 


S2, 


10 = 


= {r2,r6] 


■ (^^4,^8) 


■ (rio,ri4) ■ (ri2,ri6) ■ S2, 




11 = 


= ir2,rQ] 




■ (^"9,^3) ■ (^^10,^4) ■ (ri2,ri6) 


• S2. 



These constructions, as we saw before in Section 11.21 are not unique. A different 
construction of F2 is given in Figure 13.31 
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Figure 3.3: Another construction of Benjamin Franklin's 8x8 square F2. 



Interestingly, the Hilbert basis of 8 x 8 pandiagonal Franklin squares is a subset of 
the Hilbert basis of 8 x 8 Franklin squares. The thirty-two squares generated by the 
action of the group G on T2 and /2-T2 form the Hilbert basis of 8 x 8 pandiagonal 
Franklin squares. The pandiagonal Franklin squares in Figure 13.41 were constructed 
by Ray Hagstorm using the minimal Hilbert basis of Pandiagonal Franklin squares 
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Figure 3.4: Pandiagonal Franklin squares constructed by Ray Hagstorm |36|. 



Let Fs{s) denote the number of 8 x 8 Franklin squares with magic sum s. We 
used the program CoCoA to compute the Hilbert-Poincare series JZ^o^sis)^^ and 
obtained 

{(t36 _ ^34 _^ 28 + 33 130 ^ 233 + 390 + 947 + 1327^22 + 1991 ^20 
+ 1878^1^ + + 1327 + 947 1^^ ^ 399^10 + 233 + 33^^ + 28 

-t2 + l)}/{(t2_ 1)7(^6 _ 1)3(^2 ^1)6| 

= 1 + 34 + 64 + 483 + 1152 + 4228 + 9792 1^^ + 25957 1^^ + ■ ■ ■ 

We recover the Hilbert function Fgls) from the Hilbert-Poincare series by interpo- 
lation (see Section [TT^ . The formulas for the number of 8 x 8 pandiagonal Franklin 
squares in Theorem 11.4.41 are derived similarly. 

Natural 8x8 Franklin squares always have magic sum 260. From Theorem 11.4.31 
we find that -F8(260) is 228,881,701,845,346. This number is an upper bound for the 
number of natural 8x8 Franklin squares. The actual number of such squares is still 
an open question. 



3.2 A few aspects of 16 x 16 Franklin squares. 

Finding the minimal Hilbert basis for the cone of 16 x 16 Franklin squares is compu- 
tationally challenging and remains an unresolved problem. However, we can provide 
a partial Hilbert basis that enables us to construct Benjamin Franklin's square F3, as 
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well as the square N4. The following lemma proves that every 8x8 Franklin square 
corresponds to a 16 x 16 Franklin square. 

Lemma 3.2.1. Let M be an 8 x 8 Franklin square. Then the square T constructed 
using M as blocks (as in Figure fX3)) zs a 16 x 16 Franklin square. 



T = 



M 


M 


M 


M 



Figure 3.5: Constructing a 16 x 16 Franklin square T using an 8 x 8 Franklin square M. 



Proof. Let the magic sum of M be s. The half-columns and half- rows of T add up to 
s since they are the columns and rows of M, respectively. Also the columns and rows 
of T add to 2s. The bent diagonals of T sum to 2s (see Figure ITBl for an explanation). 
Since the 2x2 subsquares of M add to s/2, we infer that the 2x2 subsquares of T 
add to s/2. Thus T is a 16 x 16 Franklin square with magic sum 2s. □ 



M 



Since bend dingonal sum = magic 
A+B-tC+D+E+F+G+H = magic si 



Figure 3.6: Bent diagonals of T add to twice the magic sum of M. 



Consider the set S of 16 x 16 Franklin squares obtained by applying the symmetry 
operations listed in Theorem 11.4.11 to the squares constructed by applying Lemma 
13.2.11 to the ninety-eight elements of the minimal Hilbert basis of 8 x 8 Franklin 
squares (for example, SI in Figure 1X71 is constructed from the 8x8 Franklin square 
Tl in Figure ItTTjl and the 16 x 16 Franklin squares S2 and S3 in Figure IT71 Observe 
that one-fourth the magic sum of a 16 x 16 Franklin square is always an integer 
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because its 2 x 2 subsquares add to this number. This imphes that the squares in B 
are irreducible, for they have magic sums 8 or 12 (it is easy to verify that there are 
no 16 X 16 Frankhn squares of magic sum 4). Therefore, i3 is a subset of the minimal 
Hilbert basis for the cone of 16 x 16 Franklin squares. Thus, B forms a partial Hilbert 
basis. 



1 (I 1 1 (1 I 1 1 I (1 1 I 1 I 1) I 1 I I) 1 1 f) 1 1 1 I) i 1 I I) 1 1 

1 OlOIOIOlOIOlO 10 I OlOIOlOIOIOlO 10 lOlOlOIOlOlOIOl 
0101010101010101 10101010101010 10 10101010101010 10 
1 0101010101010 10 1 0101010101010 10 101010101010101 
OIOIOIOIOIOIOIOI OlOIOIOlOIOIOIOI OlOIOIOIOIOIOIOI 
I OlOIOIOlOIOlO 10 lOIOIOlOIOIOlOl lOIOIOlOIOIOlOl 

lOIOIOlOIOIOlOl lOIOIOlOIOIOlOl 1 OlOIOIOlOIOlO 10 

1 OlOIOIOlOIOlO 10 lOIOIOlOIOIOlOl I OlOIOIOlOIOlO 10 

lOIOIOlOIOIOlOl 1 OlOIOIOlOIOlO 10 1 OlOIOIOlOIOlO 10 

1 OlOIOIOlOIOlO 10 1 OlOIOIOlOIOlO 10 lOIOIOlOIOIOlOl 

lOIOIOlOIOIOlOl 1 OlOIOIOlOIOlO 10 I OlOIOIOlOIOlO 10 

1 OlOIOIOlOIOlO 10 1 OlOIOIOlOIOlO 10 lOIOIOlOIOIOlOl 
OIOIOIOIOIOIOIOI OIOIOIOIOIOIOIOI OIOIOIOIOIOIOIOI 
I OlOIOIOlOIOlO 10 lOIOIOlOIOIOlOl lOIOIOlOIOIOlOl 
OIOIOIOIOIOIOIOI OIOIOIOIOIOIOIOI I OlOIOIOlOIOlO 10 
1 OlOIOIOlOIOlO 10 OIOIOIOIOIOIOIOI I OlOIOIOlOIOlO 10 

SI S2 S3 



Figure 3.7: Elements of a partial Hilbert basis of 16 x 16 Franklin squares. 



We obtain F3 and N4 as follows: 

F3 = gl + 17 ■ g2 + 32 ■ g3 + g4 + 64 ■ g5 + 128 ■ g6 + 2 ■ g7 + 2 • g8 + 7 ■ g9 
+gl0 + 2-gll; 

N4 = gl + 17 ■ g2 + 64 ■ g3 + g4 + 32 ■ g5 + 128 ■ g6 + 2 ■ g7 + 2 ■ g8 + 7 ■ g9 
+gl0 + 2-gll. 

Since every 8x8 Franklin square corresponds to a 16 x 16 Franklin square by 
Lemma [3. 2. 11 the formulas for the number of 8 x 8 Franklin squares of magic sum s in 
Theorem II .4.31 also yields a lower bound for the number of 16 x 16 Franklin squares 
of magic sum 2s. 



3.3 Symmetries of Franklin Squares. 

In this section we prove Theorem 11.4.21 which asserts that the new Franklin squares 
Nl, N2, N3, and N4 are not derived from Benjamin Franklin's squares Fl, F2, or F3 
by symmetry operations. We first prove Theorem 11.4. 11 
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Rotation, reflection, and taking the transpose are plainly symmetry operations on 
Franklin squares. The proof of Theorem 11.4.11 follows from Lemmas 13. 3. H 1^3. 2t 13.3.3} 
andimi 

Let Sn denote the group of x n permutation matrices acting on n x matrices. 
As earlier, let {ri,rj) denote the operation of exchanging rows i and j of a square 
matrix, and let (cj, Cj) denote the analogous operation on columns. 

Lemma 3.3.1. Let G be the subgroup of Ss generated by 

{(ci, cs), (C5, C7), (C2, C4), (C6, Cg), (ri, rs), (rg, Vr), (rs, r^), (rg, rg)}, 
and let H be the subgroup of Sie generated by 

{(Ci, C3), (C2, C4), (C3, C5), (C4, Ce), (C5, C7)(c6, Cg), (cg, Cn), (cio, C12), 

(cii,ci3), (ci2,ci4), (ci3,ci5)(ci4,ci6), (ri,r3), (r2,r4), (r3,r5), (r4,r6), 

('^5,'^7)('^6,'^8), (^^9,^1), (^^10,^2), (^1,^3), (^2,^^14), (^3,^^15), (^4,^6)}- 

The row and column permutations from the group G map 8x8 Franklin squares to 
8x8 Franklin squares, while the row and column permutations from the group H map 
16 X 16 Franklin squares to 16 x 16 Franklin squares. 

Proof. Clearly exchanging rows or columns of a Franklin square preserves row and 
column sums. Half-row and half-column sums are preserved because the permutations 
of rows and columns included here operate in some half of a Franklin square. That 
2x2 subsquare sums are preserved follows from the fact that every alternate pair 
of entries in a pair of columns or rows add to the same sum (see Figure 13.81 for 
an explanation). For any 3x3 subsquare of a Franklin square, the two sums of 
diagonally opposite elements are equal (see Figure ITHl for details). This implies that, 
if we permute alternate rows or alternate columns then the new entries preserve bent 
diagonal sums. Observe that for the preservation of bent diagonal sums, it is critical 
that the alternate row and column permutations be restricted to act in one half of a 
Franklin square (see Figure IT. 211 for examples). □ 
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Because of the 2x2 subsquare sum property. 

we have a+b+d+e = d+e+g+h. 

Therefore a+b = g+h. 

Similarly, we get b+c = h+i. ^^^^ 

So, we get a+b+h+i = b+c+g+h. 

Therefore a+i = c+g. 

Similarly, a+d = c+f. 



When we interchange these two alternate columns 
, a gets replaced by c, and i gets replaced by g 
in the bent diagonal. 




Figure 3.8: Properties of Franklin squares. 



Lemma 3.3.2. Let S be the subgroup of Siq generated by the set 

{(ri, rg), (r2, rg), (rg, r^), {r^, rg), (rg, rig), (rio, r^), (rn, rig), {r^, rig), 

(Ci, C5), (C2, Ce), (C3, C7), (C4, Cs), (Cg, C13), (cio, C14), (cn, C15), (Ci2, Cie)}. 

r/ie roti; anc? column permutations from the subgroup S map 16 x 16 Franklin squares 
to 16 X 16 Franklin squares. 

Proof. Half-column and half-row sums are preserved because the specified row or 
column exchanges only affect some half of the 16 x 16 Franklin square. Because 
of the 2x2 subsquare sum property, we see that every 4x4 subsquare adds to 
the common magic sum. Hence the two sums of diagonally opposite elements in a 
5x5 subsquare are equal (see Figure EiH for an explanation). This implies that the 
bent diagonal sums are preserved, again because they operate in only one-half of a 
Franklin square. It is easy to verify that all other sums are preserved under the action 
of elements of 5*. □ 

Because of the 4x4 subsquare property, 
we have a+b+c+d = u+v+w+x 
and b+c+d+e = v+w+x+y. 
Therefore a+y = u+e. 



Figure 3.9: Properties of 16 x 16 Franklin squares. 



a b 
t S 



Observe that the groups G, H, and S are commutative, for each nonidentity ele- 
ment in each of these groups has order 2. Therefore the order of G is 2^, the order of 
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H is 2^4, and the order of 5" is 2^^ 

Lemma 3.3.3. The operation of interchanging the first n/2 columns (respectively, 
rows) and the last n/2 columns (respectively, rows) of an n x n Franklin square is a 
symmetry operation. 

Proof. These operations preserve half-row and half-column sums. The row sums and 
column sums do not change. Bent diagonal and 2x2 subsquare sums are preserved 
because of continuity (see Figure EUni for examples). □ 
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To get this square from Fl 
interchange the first 
four rows with the 
last four rows. 



To get this square from F3 
interchange the first 
eight rows with the 
last eight rows 
and then interchange 
the first eight columns 
with the 

last eight columns. 
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Figure 3.10: Constructing Franklin squares by simultaneous row (or column) exchanges of Franklin 
squares described in Lemma 13.3.31 



Lemma 3.3.4. Simultaneously interchanging all the adjacent columns (respectively, 
rows) i and i + 1 (i = 1,3,5, ... ,n — 1) of an n x n Franklin square is a symmetry 
operation. 

Proof. It is clear that row, column, half-row, and half-column sums are preserved 
by these operations. Moreover 2x2 subsquare sums are preserved because every 
alternate pair of entries in a pair of columns or rows add to the same sum (see Figure 
13. 8p . Bent diagonal sums are preserved because of the 2x2 subsquare sum property. 
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See Figure IH. 1 II for an explanation (the explanation for bent diagonal sums of 16 x 16 
Franklin squares is similar). □ 
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A + B + C + D + E + F + G + H = magic sum. 
From the 2x2 subsquare sum property we get: 
A+B+C+D+E+F+G+H+ 
a + b + c + d + e + f+ g + h = twice magic sum. 
Therefore 

a + b + c + d + e + f+ g + h = magic sum. 



Figure 3.11: Properties of 8 x 8 Franklin squares. 

Are the symmetry operations given in Theorem 1 1 . 4 . 1 1 all the symmetry operations 
of a Franklin square? We do not know the answer to this question, but the symmetries 
described in this section enable us to interchange any two rows or any two columns 
within any half of a Franklin square and get a Franklin square. Note that certain row 
or column exchanges are not symmetry operations unless they are accompanied by 
other simultaneous row or column exchanges. We now show that Nl, N2, N3, and 
N4 are not symmetric transformations of Fl, F2, or F3. 

Lemma 3.3.5. The squares Fl and F2 can he transformed by means of symmetry 
operations neither to each other nor to any of the nonisomorphic squares Nl, N2, or 
N3. 

Proof. By definition, symmetry operations map a Franklin square to another Franklin 
square. We can permute the entries of the Franklin square F2 to get Fl, Nl, N2, 
and N3 (see Figure I3.12|) . The permutation that maps F2 to N2 is not, however, a 
symmetry operation: Fl does not transform to a Franklin square under this permu- 
tation since bent diagonal sums are not preserved. The other permutations of F2 in 
Figure 13.121 likewise fail to be symmetries. Again, Fl does not map to a Franklin 
square under these permutations for half-column sums are not preserved. Thus F2 
cannot be transformed to Fl, Nl, N2, or N3 using symmetry operations. Similarly, 
the permutations of the entries of Fl that map it to F2, Nl, N2, and N3, respec- 
tively, are not symmetry operations because F2 is not mapped to a Franklin square 
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under any of these permutations (in these instances half- row sums are not preserved). 
The permutations that map the square Nl to N2 and N3, and the permutations that 
map N2 to Nl and N3 are not symmetry operations because Fl is not mapped to a 
Frankhn square under these permutations (half-column sums of Fl are not preserved 
for all these permutations). Similarly, the permutations that map N3 to Nl and N2 
are not symmetry operations because F2 is not mapped to a Franklin square under 
these operations. Therefore, the squares Nl, N2, and N3 are not isomorphic to each 
other. □ 



Fl F2 Nl 



g4 


h4 


e4 


14 


c4 


d4 


a4 


b4 




al 


a2 


a3 


a4 


a5 


a6 


a7 


a8 




h3 


h4 


hi 


h2 


h7 


h8 


h5 


h6 


S3 


h3 


e3 


B 


c3 


d3 


a3 


b3 


bl 


b2 


b3 


b4 


b5 


b6 


b7 


b8 


bl 


b2 


b3 


b4 


b5 


b6 


b7 


bS 


e5 


f5 


g5 


h5 




b5 


c5 


d5 


cl 


c2 


c3 


c4 


c5 


c6 


c7 


c8 


cl 


c2 


c3 


c4 


c5 


c6 


c7 


c8 


e6 


f6 


ofl 


li6 


a6 


b6 


c6 


d6 


dl 


d2 


d3 


d4 


d5 


d6 


d7 


dS 


e3 


e4 


el 


e2 


e7 


eS 


e5 


e6 


f7 


el 


h7 


g7 


b7 


ii7 


d7 


c7 




el 


e2 


e3 


e4 


e5 


e6 


e7 


e8 




d3 


d4 


dl 


d2 


d7 


d8 


d5 


d6 


f8 


eS 


hS 


gS 


b8 


ii8 


dS 


c8 


fl 


f2 


13 


14 


f5 


f6 


17 


IS 


fl 


n 


13 


14 


15 


16 


17 


IS 


h2 


§2 


f2 


e2 


d2 


c2 


b2 


a2 


gl 


g2 


Ll3 


g4 


g5 


g6 


g7 


g8 


gl 


g2 


g3 


g4 


g5 


g6 


g7 


gS 


hi 


gl 


fl 


el 


dl 


cl 


bl 


al 


hi 


h2 


h3 


h4 


h5 


h6 


h7 


h8 


a3 


a4 


al 




a7 


a8 


a5 


a6 






N2 












N3 








al 


a2 


a3 


a4 


dS 


d7 


d6 


d5 




d5 


h4 


e4 


a5 


h5 


d4 


a4 


e5 




bl 


b2 


b3 


b4 


c8 


c7 


c6 


c5 


d6 


h3 


e3 


a6 


h6 


d3 


a3 


e6 


cl 


c2 


c3 


c4 


bS 


b7 


b6 


b5 


b4 


15 


g5 


c4 


f4 


b5 


c5 


g4 


dl 


d2 


d3 


d4 


aS 


a7 


a6 


a5 


b3 


16 


g6 


c3 


B 


b6 


c6 


g3 


el 


e2 


e3 


e4 


hS 


h7 


h6 


h5 




e7 


h7 


d2 


e2 


a7 


d7 


h2 


11 


12 


B 


14 


g8 


g7 


g6 


g5 


al 


e8 


hS 


dl 


el 


aS 


d8 


hi 


gl 


g2 


g3 


g4 


18 


17 


16 


13 


c7 




12 


b7 


g7 




b2 


17 


hi 


h2 


h3 


h4 


eS 


e7 


e6 


e5 


c8 


gl 


fl 


b8 


gS 


cl 


bl 


18 



Figure 3.12: The abstract permutation of F2 that gives Fl, Nl, N2, and N3. 



Lemma 3.3.6. The square F3 cannot he transformed to N4 using symmetry opera- 
tions. 

Proof. Permuting the entries of F3 to get N4 is achieved by simultaneously inter- 
changing columns 1 and 15, columns 2 and 16, columns 7 and 9, and columns 8 and 
10 of F3. This permutation is not a symmetry operation. To see this, note that 
the square A obtained by transposing F3 is a Franklin square f Theorem 11.4. 1|) . But 
(ci, Ci5)(c2, Ci6)(c7, C9)(c8, Cio)-A is uot a Franklin square, since bent diagonal sums 
are not preserved. □ 

Lemmas 13 . 3 . 51 and 13 . 3 . 6t in tandem, establish Theorem 11.4.21 
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Chapter 4 

Symmetric Magic Squares and the 
Magic Graphs Connection 

What immortal hand or eye 
Dare frame thy fearful symmetry? 
- William Blake. 

4.1 Hilbert bases of polyhedral cones of magic labelings. 

In this section we derive some results about Hilbert bases of cones of magic labelings 
of graphs and also prove Proposition II. 5. HI 

Lemma 4.1.1. Let G be a graph with n vertices. A labeling L of G with magic sum 
s can be lifted to a magic labeling L' ofVn with magic sum s. 

Proof. Since G is a subgraph of r„, every labeling L of G can be lifted to a labeling 
L' of r„, where 

{L{eij) if Cij is also an edge of G, 
otherwise. 

Since the edges with nonzero labels are the same for both L and L', it follows that 
the magic sums are also the same. □ 

Lemma 4.1.2. Let G be a graph with n vertices. The minimal Hilbert basis of Gq 
can be lifted to a subset of the minimal Hilbert basis ofG^^. 
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Proof. If L is an irreducible magic labeling of G, then clearly it lifts to an irreducible 
magic labeling L' of r„. Since the minimal Hilbert basis is the set of all irreducible 
magic labelings, we get that the minimal Hilbert basis of Cq corresponds to a subset 
of the minimal Hilbert basis of Cr . □ 

-L n 

For example, the magic labelings 01 and 02 of the octahedral graph in Figure ^31 
correspond to the magic labelings al and a2, respectively, of Fg (see Figure WH^ . 
Similarly, we can prove: 

Lemma 4.1.3. For a digraph D with n vertices, a magic labeling L with magic sum 
s can be lifted to a magic labeling L' of Tin with the same magic sum s. The minimal 
Hilbert basis of Cu can be lifted to a subset of the minimal Hilbert basis ofCun- 

Lemma 4.1.4. Let D be a digraph with n vertices. All the elements of the minimal 
Hilbert basis of Cu have magic sum 1. 

Proof. It is well-known that the minimal Hilbert basis of semi-magic squares are 
the permutation matrices (see jHI]) and therefore have magic sum 1. The one-to-one 
correspondence between magic labelings of H„ and semi-magic squares implies that 
the minimal Hilbert basis elements of Cn„ have magic sum 1. It follows by Lemma 
14. 1.31 that all the elements of the minimal Hilbert basis oi Cd have magic sum 1. □ 

We now present the proof of the fact that perfect matchings of bipartite graphs 
correspond to the elements of the minimal Hilbert basis of its corresponding digraph. 
Proof of Proposition \1.5.c\ 

Every element of the minimal Hilbert basis of corresponds to a perfect match- 
ing of B by Lemma 14.1.41 Moreover, all the magic labelings of of magic sum 1 
belong to the minimal Hilbert basis of because they are irreducible. Since per- 
fect matchings of B are in one-to-one correspondence with magic labelings of B^ of 
magic sum 1, we derive that perfect matchings of B are in one-to-one correspondence 
with the elements of the minimal Hilbert basis of Cb^- It follows that Hbj^{1) is the 
number of perfect matchings of -B. □ 

For example, consider the Octahedral graph with the given orientation Dq in 
Figure 11.271 The minimal Hilbert basis of Dq is given in Figure 14.11 The perfect 
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matchings of the bipartite graph Goq corresponding to the minimal Hilbert basis 
elements of Co^ is given in Figure We derive Hu^{r) = r + 1 and thereby verify 
that the number of perfect matchings of Gn^ is indeed 2. 




Figure 4.1: The minimal Hilbert basis of the cone of magic labelings of Do- 




Figure 4.2: Perfect matchings of Gdq corresponding to the minimal Hilbert basis elements of Dq 
(see Figure 



4.2 Counting isomorphic simple labelings and Invariant rings. 

Let Sn denote the group of permutations that acts on the vertex set {vi,V2, ■ ■ ■ ,Vn} 
of G. Let Cij denote an edge between the vertices Vi and vj. The action of Sn on the 
vertices of G translates to an action on the labels of the edges of G by 

(T{L{eij)) = L{e^(^iy(j)) where a e S„. 

Two labelings L and L' of G are isomorphic if there exists a permutation cr in Sn 
such that L'{eij) = L{e„(^i)a-{j)) , i-e. cr(L) = L'. A set S = {gi,g2, ■ ■ ■ ,gr} is said to 
generate the minimal Hilbert basis of the cone of magic labelings of G if every element 
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Figure 4.3: The minimal Hilbert basis of magic labelings of the tetrahedral Graph. 

of the minimal Hilbert basis is isomorphic to some gi in 5*. For example, Tl generates 
the minimal Hilbert basis of the cone of magic labelings of the tetrahedral graph (see 
Figure ESI)- Observe that we get T2 by permuting the vertices V2 and of Tl, and 
T3 by permuting the vertices W3 and of Tl. 

Proposition 4.2.1. Let L he a magic labeling in the minimal Hilbert basis of the 
cone of magic labelings of a graph (or a digraph), then all the labelings isomorphic to 
L also belong to the minimal Hilbert basis. 

Proof. Let L' be a labeling isomorphic to L and a in Sn be such that (j{L) = L' . 
Suppose V does not belong to the Hilbert basis. Then V is reducible and can be 
written as sum of two labelings: L' = Li + L2. But, a^^{L') = L. Therefore, 
(T^^{Li) + (T~^(L2) = L. This is not possible because L is irreducible, since it belongs 
to the minimal Hilbert basis. Therefore, we conclude that L' must also belong to the 
Hilbert basis. □ 

A labeling of G is called a simple labeling if the labels are or 1. Invariant theory 
[HU] provides an efficient algebraic method of counting isomorphic simple labelings of 
a graph G. Let L be a simple labeling of G. Let denote the monomial 

i,j=l,...,n 

Consider the polynomial 

(X^)®= ^X'^(^), where X'^(^)= J] 

a&S„ i,j=l,...,n 




CHAPTER 4. The Magic Squares and Magic Graphs Connection 60 

Observe that (X^) is an invariant polynomial under the action of Sn on the 
indices of the variables Xij. Let k be any field. The set of polynomials invariant in 
the polynomial ring k[xij] under the action of the group Sn is called the invariant 
ring of 5*^ and is denoted by k[xij]'^". See [59J for an introduction to invariant rings. 

Consider the simple labeling Lq of r„ associated to G: 

{1 if Cij is an edge of G, 
otherwise. 

Then the polynomial (X^)® evaluated at Lq counts the number of labelings of G 
that are isomorphic to L. 



[8] 2 1 [4] 2 




Edge Label 



01 02 



Figure 4.4: Generators of the minimal Hilbert basis of magic labelings of the Octahedral Graph. 

For example, consider the Octahedral graph O with the labeling 01 in Figure 031 
Then, 

{X'^^^ = Xi4X23X;iQ + X1QX23X45 + 0:13X263^45 + ^12X3^X45 + X15X23X4Q + 0:13X25X46 + 
3^162^252^34 + 2:15X262:34 + X12X34X56 + X12X35X46 + Xi3X242:56 + 2:16X242:35 + 
2^142:262^35 + 2:15X242:36 + 2:14X25X36. 

Substituting X12 = X13 = X14 = X15 = X23 = 2:24 = 2:26 = 2:35 = X36 = x^4.5 = = 
X56 = 1 and xii = X22 = 2:33 = X44 = 2:55 = X66 = 2:16 = X25 = X34 = in (X*^^)®, we 
get that 

Therefore, there are 8 magic labelings in the 5*6 orbit of the magic labeling 01 of 
the octahedral graph. Similarly, there are four magic labelings in the orbit of 02. 
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The generators of the Hilbert basis of the Octahedral graph are given in Figure 14.41 
The numbers in square brackets in the figures indicate the number of elements in the 
orbit class of each generator throughout the article. 

For digraphs, we assign a variable Xij to every directed edge etj, and use the 
corresponding invariant ring to count isomorphic simple labelings. See jHI] for more 
aspects of labeled graph isomorphisms and invariant rings. 

Since all the elements of the minimal Hilbert basis of Cd, where D is a digraph, 
are simple labelings, we can use invariant theory effectively. Since the number of 
elements in the minimal Hilbert basis of semi-magic squares (and hence Cn„) is n! 
(see [51]), we can list all the generators of the minimal Hilbert basis of Cu„- The 
generators of the Hilbert basis of Cng are given in Figure 14.51 
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Figure 4.5: Generators of the minimal Hilbert basis of 6 x 6 semi-magic squares. 



CHAPTER 4. The Magic Squares and Magic Graphs Connection 



62 



Algorithm 4.2.1. Computing minimal Hilbert basis of a finite digraph D with n 
vertices. 

Input: A digraph D with n vertices and the set of n x n permutation matrices. 

Output: The minimal Hilbert basis of the finite digraph D. 

Step 0. List a set of generators of the minimal Hilbert basis of Cn„. 

Step 1. Choose all the elements hi among the generators of the minimal Hilbert 
basis of Cn„ which have the edges not in D labeled 0. Delete the edges in hi that are 
not in D to get a magic labeling Qi of D. 

Step 2. Qi and the magic labelings isomorphic to gi form the minimal Hilbert basis 
of the cone of magic labelings of D. 

For example, consider the digraph D that have all the edges of Hg except the 
loops. Then the minimal Hilbert basis of D are the 265 labelings corresponding to 
the labelings si, s2, s7, and sll in Figure and their isomorphic magic labelings. 

4.3 Polytopes of magic labelings. 

The proofs of Theorems 11.5.41 ll.5.5t and ll.5.5l and Corollarv 11.5.4.11 are presented in 
this section. 

Let G be a positive graph. An element (3 in the semigroup Scq is said to be 
completely fundamental, if for any positive integer n and a, a' G 5'cg, njS = a + a' 
implies a = ij3 and a' = {n — i)(3, for some positive integer i, such that < i < n 
(see ISl). 

Lemma 4.3.1. Vg is a rational polytope. 

Proof. Proposition 4.6.10 of Chapter 4 in jHlj states that the set of extreme rays 
of a cone and the set of completely fundamental solutions are identical. Proposition 
2.7 in [S3 states that every completely fundamental magic labeling of a graph G has 
magic sum 1 or 2. Thus, the extreme rays of the cone of magic labelings of a graph 
G are irreducible 2-matchings of G. We get a vertex of Vg by dividing the entries of 
a extreme ray by its magic sum. Thus, Vg is a rational polytope. □ 
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Lemma 4.3.2. The dimension ofVc is q — n + b, where q is the number of edges of 
G, n is the number of vertices, and b is the number of connected components that are 
bipartite. 

Proof. Ehrhart's theorem states that the degree of Hair) is the dimension of Vg [12]. 
The degree of Hair) is g — n + 6 by Theorem 11.5.21 Therefore, the dimension of Vg 
is q — n + b. □ 

Lemma 4.3.3. The d-dimensional faces of Vg are the d-dimensional polytopes of 
magic labelings of positive subgraphs of G with n vertices and at most n — b + d edges. 

Proof. An edge e labeled with a zero in a magic labeling L of G does not contribute 
to the magic sum, therefore, we can consider L as a magic labeling of a subgraph 
of G with the edge e deleted. Since a face of Vg is the set of magic labelings of G 
where some edges are always labeled zero, it follows that the face is also the set of 
all the magic labelings of a subgraph of G with these edges deleted. Similarly, every 
magic labeling of a subgraph H with n vertices corresponds to a magic labeling of 
G, where the missing edges of G in if are labeled with 0. Now, let if be a subgraph 
such that the edg labeled zero for every magic labeling of H. Then 

the face defined by H is same as the face defined by the positive graph we get from 
H after deleting the edges Cj-i, . . . , Crm- Therefore, the faces of Vg are polytopes of 
magic labelings of positive subgraphs. 

By Lemma r4.3.2t the dimension of Vg is q—n+b. Therefore, to get a d-dimensional 
polytope, we need to label at least q — n + b — d of G edges always 0. This implies 
that the o?-dimensional face is the set of magic labelings of a positive subgraph of G 
with n vertices and at most n — b + d edges. □ 

The proof of Theorem 11.5.41 follows from Lemmas 14.3.11 14.3.2[ and 14.3.31 We can 
now prove CoroUarv 11.5.4.11 
Proof of Corollary \1.5.4.i 

It is clear from the one-to-one correspondence between magic labelings of and 
symmetric magic squares that Vr„ has the following description: 
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Pr„ = {i^ = (i^(e.,)GM^;L(e,,)>0;l<^,J<r^,^<J, 
E5=i H^ji) + EJ=i+i = 1 for i = 1, . . . , n}. 

Since the graph r„ has "■^"'^^^ edges and n vertices, and every graph is a subgraph 
of r^, it foUows from Theorem 11.5.41 that the dimension of Vr^ is ; the d- 

dimensional faces of Pr„ are d-dimensional polytopes of magic labehngs of positive 
graphs with n vertices and at most n + d edges. 

We can partition the vertices of into two equal sets A and B in (^"^^) ways: 
Fix the vertex vi to be in the set A, then we can choose the n vertices for the set B 
in (^"^^) ways, and the remaining n — 1 vertices will belong to the set A. By adding 
the required edges, we get a complete bipartite graph for every such partition of the 
vertices of r2n- Thus, the number of subgraphs of that are isomorphic to Kn^n is 
^2n-iy Therefore, there are (^"^^) faces of 'Pr2n ^^at are Birkhoff polytopes because 
every isomorphic copy of Kn^n contributes to a face of 'Pr2„- ^ 

We now prove our results about polytopes of magic digraphs. 
Proof of Theorem ] 1.5. 51 By Lemma f4. 1.41 all the elements of the Hilbert basis of Cd 
have magic sum 1. Since the extreme rays are a subset of the Hilbert basis elements, 
it follows that the vertices of Vd are integral. Since Vd = 'Pgdi follows by Theorem 
11.5.41 that the dimension of Vd is q — 2n + b; the (i-dimensional faces of Vd are the d- 
dimensional polytopes of magic labelings of positive subdigraphs of D with n vertices 
and at most 2n — b + d edges. □ 

We derive our results about the faces of the Birkhoff polytope as a consequence. 
Proof of Theorem ] 1.5. (A The one-to-one correspondence between semi- magic squares 
and magic labelings of n„ gives us that Vn„ = Bn. Since every digraph with n vertices 
is a subdigraph of n„, by Theorem ll.5.5| it follows that its (i-dimensional faces are 
(i-dimensional polytopes of magic labelings of positive digraphs with n vertices and 
at most 2n ~ 1 + d edges. Since the vertex set of a face of Bn is a subset of the 
vertex set of i?„ it follows that the vertices of Vd, where D is a positive digraph, are 
permutation matrices. □ 
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Our results enable us to reprove some known facts about the Birkhoff polytope as 
well. For example, Theorem 11.5.51 gives us that the dimension of i?„ is {n — 1)^. The 
leading coefficient of the Ehrhart polynomial of -B„ is the volume of B„. This number 
has been computed for n = 1, 2, . . . , 9 (see yjj and 

4.4 Computational results 

We will now list our computational results. The numbers in square brackets in the 
figures represent the number of elements in the orbit class of each generator. 

4.4.1 Symmetric magic squares. 

The generators of the minimal Hilbert basis of Cr„ forn = 1, 2, 3, 4, 5, and 6 are given 
in Figures 14.61 14.71 14.81 14.91 14.1U1 and 14.111 respectively. It is interesting that all the 
elements of the minimal Hilbert basis of Cr„ for n = 1, . . . , 6 are 2-matchings. The 
minimal Hilbert basis elements are not, in general, 2-matchings for all n (see Figure 
14.121 for examples of irreducible magic labelings of magic sum 3). In fact, it follows 
from the results of chapter 11 in [Hj, that there exists a graph with an irreducible 
magic labeling of magic sum r if and only if r is 2 or r is odd. By Proposition 14. 1.^?} 
this implies, that there is a minimal Hilbert basis element of magic sum r of Cr„, for 
some n, if and only if, r is 2 or r is odd. The program 4ti2 was also able to compute 
the minimal Hilbert bases of C-py and C^g- 

Recall that the number of n x n symmetric magic squares is the same as HY^{r) 
(the generating functions for H^^lr) for n up to 5 are given in jSHl)- The volume of 
Vr„ is the leading coefficient of H-p^{r). 

m 



Figure 4.6: The minimal Hilbert basis of 1 x 1 symmetric magic squares. 
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Figure 4.7: Generators of the minimal Hilbert basis of 2 x 2 symmetric magic squares. 



[3] 



[1] 



2# • 3 2« • 3 2# • 3 




[1] 



Edge 


Label 
1 








'5 

Figure 4.8: Generators of the minimal Hilbert basis of 3 x 3 symmetric magic squares. 



Hr,{r) = 1 for all r > 0. 



Hr^{r) = r + 1 for all r > 0. 



1^,3 _^ ^ 7^ _^ jf 2 divides r, 



i^^ + 1^^ + + i otherwise. 



^r^ + ir^ + flfr^ + fr^ + f + |r + 1 if 2 divides r, 



+ + ii!^^ + + f + 3^ + 16 otherwise. 
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[1] 

« ^ 



[3] 



[6] 



[4] 




3 Edge Ubel 



Figure 4.9: Generators of the minimal Hilbert basis of 4 x 4 symmetric magic squares. 
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[10] 



0- 
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Figure 4.10: Generators of the minimal Hilbert basis of 5 x 5 symmetric magic squares. 



365 10 _|_ 9125 9 i 22553 8 i 55085 7 i 11083 6 i 7945 5 i 1978913 4 
3096576 3096576 688128 258048 12288 3072 387072 



I 335065 3 I 50329 2 , 1177 , -i 
48384 8064 336 



if 2 divides r, 



365 10 I 9125 9 i 22553 8 i 55085 7 i 11083 6 i 63545 5 i 15807679 4 
3096576 3096576 688128 258048 12288 24576 3096576 



I 5329855 3 i 6327137 2 , 1139917 , 27213 
774144 ' 1032192 ~'~ 344064 ~'~ 32768 



otherwise. 
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[15] 



[10] 



[45] 



[72] 





[1] 



[60] 



[60] 



» 

' 3 



[20] 



[15] 



[60] 




Edge Label 



Figure 4.11: Generators of the minimal Hilbert basis of 6 x 6 symmetric magic squares. 



4.4.2 Pandiagonal symmetric magic squares. 

Pandiagonal symmetric magic squares are symmetric magic squares sucli tliat all the 
pandiagonals also add to the magic sum (see Figure ll.4|l . The generators of the 
Hilbert basis of n x n pandiagonal symmetric magic squares for n = 3, 4, and 5 
are given in Figures 14.131 I4.14[ and 14.151 respectively. The Hilbert basis of 6 x 6 
pandiagonal symmetric magic squares contain 4927 elements and can be computed 
using the program 4ti2. 

Let Pn{r) denote the number of n x n pandiagonal symmetric magic squares with 
magic sum r. We derive: 
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Edge Label 

3 

2 

1 





Figure 4.12: Minimal Hilbert basis elements of Cr^ and Cpg of magic sum 3. 




Figure 4.13: Generators of the minimal Hilbert basis of 3 x 3 pandiagonal symmetric magic squares. 



1 if 3 divides r, 



otherwise. 



+ |r + 1 if 4 divides r, 



otherwise. 



+ i-t^ + §t^ + §t + l if 2 divides r, 



1 y.4 5_j.2 I _3_ 

384'' 192'' 128 



otherwise. 



4.4.3 Magic labelings of Complete Graphs. 



The minimal Hilbert basis of the cone of magic labelings of the complete graph Kn 
corresponds to the set of elements of the minimal Hilbert basis of Cu„ for which 
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[1] 



[4] 





Figure 4.14: Generators of the minimal Hilbert basis of 4 x 4 pandiagonal symmetric magic squares. 




[5] 

'9 2 



Edge Labd 



^ r • 

^6 



Figure 4.15: Generators of the minimal Hilbert basis of 5 x 5 pandiagonal symmetric magic squares. 



all the loops are labeled with a 0. Vk^ is an dimensional polytope with the 



description: 



Hk,{t) = 1 

1 if 2 divides r, 



otherwise. 



Jgr^ + + if r3 + fr^ + ||r + 1 if 2 divides r, 



Hk,{t) = { 



otherwise. 
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19 „9 _|_ 19 „8 _|_ 143 7 i 5_ 6 , 4567 5 i 785 4 i 10919 3 
120960 ~'~ 5376 ~'~ 4032 ~'~ 24 ~'~ 5760 ~'~ 384 ~'~ 3024 



I 955 ^2 I 857^ i i 
"■"224 280 ~'~ 



if 2 divides r, 



19 9 I J9_ 8 I JA3_ 7 I 5_ 6 , 4567 5 , 785 4 , 10919 3 



120960 ' 5376 ' 4032 '24 ' 5760 '384 ' 3024 



I 955 2 I 857 , 251 
"•"224 280 ~'~ 256 



otherwise. 



See jSH] for more aspects of magic labelings of complete graphs Kn- 



4.4.4 Magic labelings of the Petersen graph. 

The generators of the Hilbert basis are given in Figure 14.161 (numbers in square 
brackets are the number of elements in the orbit of the generators). 

[6] 




Figure 4.16: Generators of the minimal Hilbert basis of magic labelings of the Petersen Graph. 



Let Hpetersenif) dcuotc the uumbcr of magic labelings of the Petersen graph with 
magic sum r. The generating function F{t) for the Petersen graph (F(t) is also derived 
in |Sn|) is: 

t^ + t^ + 6t'^ + t+l 



Fit) 



(i-tni+t) 

Therefore, we get 



l + 6t+27t2 + 87t3 + 228t^ + 513t^ + 1034t6 + 1914t^ + . . . 
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25^3 _^ _^ _^ g divides r, 



Hpeterseny) ^ 



I, ^''^^ + h^'^ + H''^^ + T'^^ + fl'^ + tI otherwise. 
4.4.5 Magic labelings of the Platonic graphs. 
Magic labelings of the Tetrahedral Graph. 

The minimal Hilbert basis of the cone of magic labelings of the tetrahedral graph is 
given in Figure 14.31 Since all the elements of the minimal Hilbert basis have magic 
sum 1, it follows that the vertices of the polytope of magic labelings of the tetrahedral 
graph are integral points. Therefore, we get that Htetrahedralir) is a polynomial, where 
Htetrahedralir) dcuotcs the uumbcr of magic labelings of the Tetrahedral graph. We 
derive an explicit formula. 

1 3 

Htetrahedralir) = -r^ + -r + 1. 

Theorem 11.5.31 states that if a graph G is bipartite then Heir) is a polynomial. 
Thus, the tetrahedral graph is an example that proves that Hair) being a polynomial 
does not imply that G is bipartite. 

Magic labelings of the Cubical graph. 

The minimal Hilbert basis of the cone of the magic labelings of the cubical graph is 
the set consisting of CI in Figure 14.171 and the eight magic labelings isomorphic to 
CI. 

[9] 



Edge 


label 





1 



CI 



Figure 4.17: Generator of the minimal Hilbert basis of magic labelings of the Cubical Graph. 
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Let Hcubeij) denote the number of magic labelings of the Cubical graph of magic 
sum r. Observe that the cubical graph is bipartite, therefore Theorem 11.5.31 applies, 
and Hcuheij) is a polynomial. Hcubeij) is also derived in 



cubeK J 15 ^2 3 30 



The magic labelings of the Octahedral graph. 

There are 12 elements in the minimal Hilbert basis of the cone of magic labelings 
of the Octahedral graph: the 8 magic labelings in the 5*6 orbit of the magic labeling 
01, and the four magic labelings in the orbit of 02 (see Figure lOj) . Let Hodahedrai 
denote the number of magic labelings of the Octahedral graph of magic sum r. The 
generating function of Hoctahedrai{r) is given in 



H octahedral ('^ ) 



_i_^6 ^ ^ 11^4 ^ + 38^2 ^ ^ ^ 2 dividcs r. 



_i_^6 ^ ^ 11^4 ^ + 38^2 ^ ^ 15 otherwise. 



The magic labelings of the Dodecahedral graph. 

The generators of the minimal Hilbert basis of the Dodecahedral graph are given in 
Figure KT^ 





Figure 4.18: Generators of the minimal Hilbert basis of magic labelings of the dodecahedral Graph. 



Let H dodecahedral (r) dcuotc the uumbcr of magic labelings of the Dodecahedral 
graph with magic sum r. We derive: 
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H dodecahedral 



47 10 I _¥7_ 9 I ^25_ 8 , 9_ 7 , 3361 6 , 255 5 , 27625 4 
40320 ~'~ 2688 ~'~ 1792 16 ~'~ 1920 64 4032 



I 1513 3 _|_ 4691 2 i 9 , i 
~'~ 168 ~'~ 560 2 ~'~ 



if 2 divides r, 



47 10 _|_ _^rf,<i _|_ ^25_ 8 _|_ 9_ 7 , 3361 6 i 255 5 i 27625 4 
40320 ~'~ 2688 ~'~ 1792 ~'~ 16 ~'~ 1920 ~'~ 64 ~'~ 4032 



I 1513 3 _|_ 4691 2 , 567 , 229 
~'~ 168 ~'~ 560 ~'~ 128 ~'~ 256 



otherwise. 



The magic labelings of the Icosahedral graph. 



There are 4195 elements in the minimal Hilbert basis of the cone of magic labelings 
of the Icosahedral graph which can be computed using 4ti2. It is interesting that 
unlike the other platonic graphs, all the minimal Hilbert basis elements are not two- 
matchings (see Figure I1.19p . The formula for the number of magic labelings of the 
Icosahedral graph remains unresolved. 





Edge 


Lilbel 

2 




1 





Figure 4.19: Icosahedral graph : minimal Hilbert basis elements of magic sum 3. 
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A.l Proof of the minimal Hilbert basis Theorem. 

The cone generated by a set X of vectors is the smallest cone containing X and is 
denoted by cone X] so 

cone X = {XiXi + .... + XkXk\k > 0; Xi, . . . , G X; Ai, . . . , > 0}. 

Proof of Theorem \1.2.1\ 

Let C be a rational polyhedral cone, generated by bi,b2, ...,bk. Without loss of 
generality 61, 62, &fc are integral vectors. Let ai, 02, be all the integral vectors 
in the polytope V: 

V = {A161 + .... + AA|0 <\i<l{i = 1, .., k)} 

Then Oi, 02, at generate C as 61, 62, &fc occur among oi, 02, at and as V is 
contained in C. We will now show that ai,a2, ...,0^ also form a Hilbert basis. Let h 
be an integral vector in C. Then there are /xi, /i2, /Xfc > such that 

b = nibi + fi2b2 -\ h/ifcfefc. (A.l) 

Then 

b = [/iij bi + [fi2\ &2 H h [/ifej &fc + (/ii - [fJ'i\ )bi + (/i2 - LAi2j )&2 H h (/Ufc - [/i/cj )6fc. 

Now the vector 

b - [/iiJ6i LAifcJ&fc = (/ii - LAiiJ)&i H h (/ifc - LAtfcJ)6fc (A. 2) 
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occurs among 01,02, as the left side of the Equation IA.2I is clearly integral 
and the right side belong to V. Since also bi,b2,...,bk occur among oi, 02, o^, it 
follows that IA.l1 decomposes 6 as a nonnegative integral combination of 01,02, ...,0^. 
So Oi, 02, Of form a Hilbert basis. 

Next suppose C is pointed. Consider H the set of all irreducible integral vectors. 
Then it is clear that any Hilbert basis must contain H. So H is finite because it is 
contained in V. To see that H itself is a Hilbert basis generating C, let 6 be a vector 
such that bx > if X E C\{0} {b exists because C is pointed). Suppose not every 
integral vector in C is a nonnegative integral combination of vectors in H. Let c be 
such a vector, with be as small as possible (this exists, as c must be in the set V). As 
c is not in if, c = Ci + C2 for certain nonzero integral vectors Ci and C2 in C . Then 
bci < be and bc2 < be. Therefore Ci and C2 are nonnegative integral combinations of 
vectors in and therefore c is also. □ 

A. 2 Proof of the Hilbert-Serre Theorem. 

Let C be a class of A-modules and let if be a function on C with values in Z. The 
function H is called additive if for each short exact sequence 

^ M' ^ M ^ M" ^ 

in which all the terms belong to C, we have 

H{M') - H{M) + ii'(M") = 0. 

Proposition A. 2.1 (proposition 2.11, [§]). Let ^ Mq ^ Mi ^ > M„ 

be an exaet sequenee of A-modules in whieh all the modules Mj and the kernels of all 
the homomorphisms belong to C . Then for any additive function H on C we have 

n 

Y,{-^rH{M,) = 0. 

i=0 

Proof. The proof follows because every exact sequence can be split into short exact 
sequences (see jH], Chapter 2). □ 
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For any A-module homomorphism of M into A^, we have an an exact sequence, 

O^ker(0)^M ^ A^-^coker(0)^O, 

where ker(0)— s>M is the inclusion map and A^— s>coker(0) = A^/im(0) is the natural 
homomorphism onto the quotient module j29j . 

Theorem A. 2.1 (Theorem 11.1 jg,(Hilbert,Serre)). Let A = 0^=o^n « 
graded Noetherian ring. Let A he generated as a A^-algehra by say Xi, X2, Xg, which 
are homogeneous of degrees ki, k2, .., kg (all > 0). Let H he an additive function on 
the class of all finitely- generated AQ-modules. Let M he a finitely generated A-module. 
Then the Hilhert-Poincare series of M, HM{t) = Yl'^=o H{Mn)t'^ is a rational function 
in t of the form p{t) /Yil^^{ \ — t''^), where p{t) E Z[t]. 

Proof. Let M = M„, where M„ are the graded components of M, then M„ is 
finitely generated as a Ag-module. The proof of the theorem is by induction on s, the 
number of generators of A over Aq. Start with s = 0; this means that An = for all 
n > 0, so that A = Ao, and M is a finitely-generated Aq module, hence M„ = for 
all large n. Thus HM{t) is a polynomial in this case. 

Now suppose s > and the theorem true for s — 1. Multiplication by Xg is an 
A-module homomorphism of M„ into Mn+ks > hence it gives an exact sequence, say 

O^Kn-^Mn ^ M„+fc, ^Ln+k. -^0. (A.3) 

K = ^n^n, L = 0„i^n are both finitely generated A-modules and both are 
annihilated by Xg, hence they are Ao[xi, . . . , Xs_i]-modules. Applying H to I A. 31 we 
have 

H{Kn) - H{M,^) + i/(M„+fcJ - i^(WJ = 0; 
multiplying by and summing with respect to n we get 

(1 - t''^)H{M, t) = H{L, t) - t^^H{K, t) + git), 

where g{t) is a polynomial. Applying the inductive hypothesis the result now 
follows. □ 
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In this chapter, we provide some basic algorithms to compute Hilbert bases, Hilbert- 
Poincare series, and toric ideals. A knowledge of Grobner bases is assumed. An 
excellent introduction to Grobner bases is given in [221 • Many available computer 
algebra packages (for example Maple and CocoA) can compute Grobner bases. 

B.l Algorithms to compute Hilbert bases. 

We describe Algorithm 1.4.5 in to compute the Hilbert basis of a cone Ca = {x : 
Ax = 0,x > 0}. 

Let A be an m X n matrix. We introduce 2n + m variables ti,t2, --tm, xi, ..,Xn, 
yi,y2, ■■,yn and fix any elimination monomial order such that 

{ti,t2, ..tm} > {Xl, .., Xn} > {Vl, 1/2, ••, Vn} . 

Let Ia denote the kernel of the map 

C[xi, ...,Xn,yi,...,yn]^ C[ti, . . .,tm,t];\ . . .,t;;^\yi, . . .,yn], 

m m 

xi^yiYl ^T' , • • • , 3;„ ^ ?/„ JJ fj"' ,yi^yi,...,yn^yn. 

j=i i=i 
We can compute a Hilbert basis of Ca as follows. 

Algorithm B.1.1 (Algorithm 1.4.5, 

1. Compute the reduced Grobner basis Q with respect to < for the ideal I a. 
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2. The Hilbert basis of Ca consists of all vectors [3 such that — appears in Q. 
For example, let 

1 -1 
-2 2 

To handle computations with negative exponents we introduce a new variable t 
and consider the lexicographic ordering 



A 



t > ti > t2 > Xi > X2 > yi > 2/2- 

We compute the Grobner basis of I a = (a^i — yitft'^, X2 — 2/2^2^? ^1^2^ — 1) with 
respect to the above ordering and get: 

Ia = ( X1X2 - ym , tm - 4x1, tiX2 - tly2, tlty2 - X2, tltxi - yi, tit2t - 1) 

Therefore, the Hilbert basis is {(1, 1)}. 

See inn] for more details about this algorithm. See [38^ for more effective algorithms 
to compute the Hilbert basis. 



B.2 Algorithms to compute toric ideals. 

Computing toric ideals is the biggest challenge we face in applying the methods we 
developed in this thesis. Many algorithms to compute toric ideals exist and we present 
a few of them here. 

Let A = {ai, a2, ctn} be a subset of Z*^. Consider the map 

TT : k[x] ^ k[t^'^] (B.l) 
Xi ^ t"' (B.2) 

Recall that the kernel of vr is the toric ideal of A and we denote it by The 
most basic method to compute would be the elimination method. Though this 
method is computationally expensive and not recommended, it serves as a starting 
point. Note that every vector -u G Z" can be written uniquely a.s u = — u~ where 

and u~ are non-negative and have disjoint support. 
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Algorithm B.2.1. [Algorithm 4-5, J^J- 

1. Introduce n + d + 1 variables to, ti, .., td, xi, X2, Xn- 

2. Consider any elimination order with {ti\i = 0,...,d} > {xj]j = l,...,n}. 
Compute the reduced Grobner basis G for the ideal 



3. GnA;[x] is the reduced Grobner basis for /_4 with respect to the chosen elimination 
order. 

If the lattice points have only non-negative coordinates, the variable to is un- 
necessary and we can use the ideal (xj — : z = 1, . . . , n) in the second step of the 
Algorithm iRm 

To reduce the number of variables involved in the Grobner basis computations, it 
is better to use an algorithm that operates entirely in k[xi, . . . , We now present 
such an algorithm for homogeneous ideals. Observe that all the toric ideals we face 
in our computations in this thesis are homogeneous. 

Recall that the saturation of an ideal J denoted by (J : /°°) is defined to be 



Let ker{A) G denote the integer kernel of the dxn matrix with column vectors 
ttj. With any subset C of the lattice ker{A) we associate a subideal of 



{totit2---td ~ I5 Xit""^ — t"^"*", Xnt' 




(J : /°°) = {g e k[x] : fg e J for some r G N}. 



Jc := (X 



u 



X" -.ueC). 



We now describe another algorithm to compute the toric ideal 



Algorithm B.2.2. [Algorithm 12.3 WJ^l. 



1. Find any lattice basis L for ker{A). 



2. Let Jl := (X 



u 



X"" : M G L). 



Appendix B 



81 



3. Compute a Grobner basis of (J^ : ( )°°) which is also a Grobner basis 

of the toric ideal /_4. 

From the computational point of view, computing (J^ : {xiX2 ■ ■ is the most 

demanding step. The algorithms implemented in CoCoA try to make this step efficient 
[Hj. For example, one way to compute (J^ : {xiX2 ■ ■ would be to eliminate 

t from the ideal H := Jl + {txiX2 - ■ ■a;^ — 1) but this destroys the homogeneity 
of the ideal. It is well-known that computing with homogeneous ideals have many 
advantages. Therefore, it is better to introduce a variable u whose degree is the sum 
of the degrees of the variables Xj, i = 1, . . . , n. We then compute the Grobner basis of 
the ideal H := Jl + {xiX2 ■ ■ - Xn — u) . Then a Grobner basis for {Jl : {xiX2 ■ ■ ■ a;„)°°) 
is obtained by simply substituting u = X1X2 ■ ■ ■ Xn in the Grobner basis of H. 

Another trick to improve the efficiency of the computation of saturation ideals is 
to use the fact 

{Jl : {x,X2 ■ ■ ■ x^D = ((••• {{-h : x^) : x^) . . .) : O- 

Therefore we can compute the saturations sequentially one variable at a time. 
See ^3] for other tricks. We refer the reader to [^0] for details and proofs of the 
concepts needed to develop these algorithms and other algorithms. We now illustrate 
Algorithm IB . 2 . 21 bv applying it to an example. 

Let A = {(1,1), (2, 2), (3, 3)}. Consider the matrix whose columns are the vectors 
of ^ 

1 2 3 
1 2 3 

Then kerA = {[— 2, 1, 0], [— 3, 0, 1]}. A lattice basis of kerA can be computed 
using the software Maple, and we get a basis is 

{[-1,-1,1], [-2, 1,0]}. 

Therefore Jl = {x^ — X1X2, X2 — xl) and 

(Jl : (xiXaXg)^) = (xg - X1X2, X2 - xj, xl - X1X3) 
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which is also Ia- 

Note that many available computer algebra packages including CoCoA can com- 
pute saturation of ideals. 

B.3 Algorithms to compute Hilbert Poincare series. 

In this section, we will describe a pivot-based algorithm to compute the Hilbert 
Poincare series jT^j. Variations of this algorithm is implemented in CoCoA. Let k 
be a field and R := k[xi,X2, ■..,Xr] be a graded Noetherian ring, let Xi,X2, ...yXr 
be homogeneous of degrees ki, k2, .., kr (all > 0). Let M be a finitely generated R- 
module. Let H be an additive function on the class of i?-modules with values in Z. 
Then by the Hilbert-Serre theorem, we have 

where p{t) G Z[t]. 

Let / be an ideal of R, we will denote 

Observe that we only need to calculate the numerator < I > since the denominator 
is already known. 

Let y he a. monomial of degree (cii, dr) called the pivot. The degree of the pivot 
is d = "^l^i di. Recall the definition of ideal quotients (J : /) [HOj 

(J :/) = {(? G k[x] : fg G J}. 

Consider the following short exact sequence on graded i?-modules. 

Q^Rl{I:y)^R/I^R/{I,y)^Q 
which yields (since H is additive) 

Hn/lit) = Hjl/^J^y){t) +t'^{Hjl/(^J.y)){t). 
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This implies 

< / >=< I,y> +t'^ <I ■.y> . (B.3) 
When J is a homogeneous ideal, 

HR/iif) = i^^/in(/)W> 



where in(J) denotes the ideal of initial terms of / 

The pivot y is usually chosen to be a monomial that divides a generator of / so 
that the total degrees of (/, y) and (/ : y) are lower than the total degree of /. The 
computation proceeds inductively. 

We illustrate this algorithm with an example. Let R = k[xi,X2, ■ ■ ■ ,a;„] be the 
polynomial ring. Let R = 0^gf^-Rd where each R^ is minimally generated as a k- 
vector by all the ("^^~^) monomials of degree d. Therefore, 



ri—n ri—n \ / 



r = 1/(1 -t)". 

\ U / 

d=0 d=0 ^ ' 

Therefore we get < >= 1. We will use this information to compute HB,/{i){t), 
where / = (xi, X2, . . . , Xn)- 

Let J = (x2, . . . , Xn)- Then, ( J : Xi) = J. Therefore by Equation IB. 3[ we get 

<(J,Xi) >= (l-t^^S"i)< J>. 



That is 



< Xi, X2, . . . , X„ >= (1 - t'^^S^i) < X2, . . . , X„ > . 

Now, choosing the pivot X2, X3, . . . , x„ subsequently we get 

<Xi,X2,...,X„>= J] (1 - t^eg^O < > . 
i=l,...,n 

Now since < >= 1, we get < xi, X2, . . . , x^ >= ni=i n(l ~ t'^^^^'). 
Therefore i/ij/(xi,x2,...,x„)(^) = 1- 

See ^S] for the effects of choosing different pivots in the algorithm and also for 
other algorithms. LattE uses a different algorithm from the CoCoA algorithms j2Z] ■ 
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And no grown-up will ever understand that this is a matter of so much importance! - Antoine 
de Saint-Exupery. 



C.l Constructing natural magic squares. 
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Figure C.l: Constructing natural magic squares with an odd number of cells [J]. 



When the entries of an n x n magic square are 1, 2, 3, n^, the magic square 
is called a pure magic square or a natural magic square. Methods for constructing 
natural magic squares are known for every order. One of my favorite methods of 
constructing a natural magic square with an odd number of cells is as follows: the 
numbers 1 to are written consecutively in diagonal columns as shown in Figure 
Ky.ll The numbers which are outside the center square are then transferred to the 
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empty cells on the opposite sides of the latter without changing their order to get a 
magic square. This method is said to have been originated by Bachet de Meziriac 
(see 0). 

A pair of numbers in the set {1, 2, n^} which add to ra^ + 1 is called complemen- 
tary. A method for constructing magic squares with an even number of cells using 
complementary pairs of numbers is as follows: write the number 1 . . . consecutively 
across rows. All entries except the numbers in the two main diagonals are then re- 
placed by their complements to get a magic square (see Figure in?2|) . See [7j for details 
of these methods and other methods of constructing natural magic squares. For more 
recent developments in the construction of natural magic squares, see |13j . 
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Figure C.2: Constructing magic squares with an even number of cells 0. 

C.2 Other magic figures. 

A composite magic square is a magic square composed of a series of small magic 
squares and an example is given in Table ICll A concentric magic square is a magic 
square that remains a magic square with borders removed. An example of a modifi- 
cation of a concentric magic square devised by Frierson is shown in figure IC.41 

A set of n magic circles is a numbering of the intersections of the n circles such 
that the sum over all intersections is the same constant for all circles. Consider the 
example of a die. It is commonly known that the opposite faces of a die contain 
complementary numbers that always add up to 7. Consequently any band of four 
numbers encircling a die gives a summation of 14 (see Figure RTSl A) . These bands 
form magic circles (see Figures FC.SI B or IC.5I C). A magic sphere is a sphere that 
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Table C.l: A composite magic square 0. 



42 41 40 



35 13 14 32 31 



28 21 26 16 



23 25 27 33 



24 29 22 30 



12 19 37 36 18 15 



' sum of the 3x3 square is 75 
- sum of the 5x5 square is 125 
^■sum of the 7x7 square is 175 



10 49 48 47 



Figure C.3: A concentric magic square 0. 



contains magic circles. The sphere in Figure IC.5I A is a magic sphere. 

A magic triangle is composed of three magic squares A, B, and C such that the 
square of any cell in C is equal to the sum of the squares of the corresponding cells in 
A and B. In other words the corresponding entries of the magic squares always form 
a Pythagorean triple (c^ = + 6^). C is called the hypotenuse, and A and B are 
called the legs of the magic triangle. An example is given in Figure IC.61 

A magic star is a numbering of the intersections of a set of lines that form a star 
such that the sum over every intersection is the same for each line (see Figure IHTTI for 
an example). 

A magic carpet is a magic square in which a limited range of digits is used several 
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31 81 11 



34 40 60 43 28 



42 22 54 39 



33 53 15 68 16 44 



49 29 67 14 66 24 38 



37 36 30 35 



74 45 46 47 52 



. sum of the 3x3 square is 123 



■ sum of the 5x5 square is 205 



- sum of the 7x7 square is 287 



■ sum of the 9x9 square is 369 



Figure C.4: Variation of the concentric magic square [Jj. 






A 



Figure C.5: Magic circles and magic sphere of dice [7]. 



times. 

A magic rectangle is an m x n matrix such that all its rows add to a prescribed 
common sum, and all its columns add to the another prescribed sum. Therefore, a 
magic rectangle has two magic sums: A magic column sum and a magic row sum. 
Figure IC.8I shows embeddings of magic squares in magic rectangles and vice versa 
and such patterns are called ornate magic squares {7\). 

See |7] for more examples of magic figures. Needless to say, we sure can construct 
and enumerate most of these magic figures with our methods. 



BIBLIOGRAPHY 




C 



Figure C.6: A Magic Triangle. 




Figure C.7: A Magic Star. 
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Table C.2: A magic carpet. 
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Table C.3: A magic carpet made from the magic carpet in Table lU^ 



10 13 11 14 1 



3 11 10 1 



Figure C.8: Ornate magic squares [Jj- 
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